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ABSTRACT 


A major difficulty in the practical application of linear-quadratic 
regulator theory is how to choose the weighting matrices in quadratic 
cost functions. In this work, the control system design with optimal 
weighting matrices is applied to a helicopter in the hover and approach 
phase. The weighting matrices are calculated to extremize the closed 
loop "Total System Damping" subject to constraints on the determinants. 

The extremization is really a minimization of the effects of disturbances, 
and interpreted as a compromise between the generalized system accuracy 
and the generalized system response speed. The trade-off between the 
accuracy and the response speed is adjusted by a single parameter — the 
ratio of determinants, which is the only arbitrariness left for designer's 

V 

choice. 

By this approach an objective measure can be obtained for the design 
of a control system. The measure is to be determined by the system re- 
quirements . 
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CHAPTER I 


INTRODUCTION 

This study will explore the use of linear optimal control methods 
to design the automatic controller of a hovering helicopter. Generally 
the pilots workload is rather high for stabilizing a hovering helicopter 
[Fig. 2] and keeping position in the presence of gust over a designated 
point. Furthermore, to keep an accurate position may become more diffi- 
cult if the motions [Fig, 3] among each axis are highly coupled. 

In design of the multivariable and multi input systems, the conven- 
tional techniques such as root locus, frequency response, and describing 
function analysis which are based on the single input-single output systems, 
have been employed in trial and error fashion. Design techniques using the 
optimal control theory applied to multivariable systems in the state space 
have become popular and have been applied for the flight controller design 
or optimal control of aircraft [Ref, 1§2] . 

The optimal control design technique utilizes the minimization of a 
performance criteria; one of which is a quadratic cost function, such as 
Eq* [3-1) . Adequate selection of weightings in the quadratic cost function 
yields a stable closed loop system [Ref. 3] . However, the following questions 
remain to be answered: 

How should the weighting matrices be .selected. 

IVhere are the poles of the closed loop system in s plane located when 
control technique is used. 

tVhat are the system responses for a command input. 

How is the system influenced by random gust inputs. 
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This work attempts to provide an answer to these questions and exa- 
mines two methods for selecting weighting matrices, namely; the method 
of authorities [Ref. 6] and the optimal weighting selection method (OWEM) 
[Ref. 7]. The approach used in this investigation is as follows: First, 
the equation of motion of a helicopter comparable to a CH54B are decoupled, 
and the results are confirmed by an analogue computer study. Secondly, 
for each set of decoupled equations are derived the optimal controls with 
its feedback gains. For a range of selected factor weightings, analysis 
has been made for the root loci of the closed loop system, optimal feed- 
back gains, the integral squared position error (ISE), attitude errors (IS0, 
IS4>) and controls (ISU^ to a command input. Also are investigated the root 
mean squared position errors, attitude errors and control motion for a ran- 
dom gust input. For the chosen pole location in s plane, the closed loop 
gains are determined, and command inputs response [ISE, IS6, etc.) and gust 
responses (rms errors) are calculated. 

In Chapter II, is reviewed the application of optimal linear control 
theory, and an outline of the performed work is given. Alsb" is discussed 
the areas that need further study and investigation. ' 

In Chapter III, the details of optimal controller design are described 
for the vertical and yawing motion, followed by the longitudinal and lateral 
motions of the hover helicopter. 

In Chapter IV, conclusions and recommendations on the optimal flight con- 
troller design using these methods are discussed. 
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In the appendices are reviewed some of the theoretical aspects of 
linear optimal control theory and model following. It appears that the 
OWEM design method is also applicable to model following, as outlined in 
Chapter II. Theoretical work has been extracted from Ref. [7], applications 
and calculations from Ref . [34] and ^h^Appendix [6] from Ref. [35]. 
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CHAPTER 2 


REVIEW OF THE APPLICATION OF LINEAR OPTIMAL CONTROL THEORY 

The purpose of this study is to apply optimal control methods for a 
controller design for helicopter at hover and in the approach phase. The 
design procedures in this report utilize first the decoupling of the heli- 
copter dynamics (Appendix B) . The system to be controlled consists of the 
decoupled helicopter dynamic equations. 

In state variable form, the equations of motion are given by: 

x(t) = Ax(t) + Bu(t) (2-1) 

In most cases, x (n dimensional state vector) , represents the state 
variables of motion of a hovering helicopter, which are, respectively, x, 
u, 8, q, z, w, y, V, <J), p, and r. The control variable u (m dimensional 
vector), represents control inputs which are respectively B^^^, A^^, 6^, and 
6^. In a few cases, the state vector x also can include the output of servo 
controller and so on. 

The output vector y(t) of helicopter is assumed to be the state vector 
x(t) itself. 


y(t) = Hx(t) = x(t) 

H = I : n X n identity matrix 
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The advantages of state variable formulation is significant in 
the case of multi-variable and multi -input systems. 

The design nethod uses an optimal control method with a quadratic 
performance index that minimizes the performance index given by the 
following equation: 

J = j f Qx(t) + u'^Ct) Ru(t)) dt C2-3) 


The matrices Q and R in Eq. (2-3) are respectively n by n and n by 
m weighting matrices which place relative weights on quadratic quantities 
x^, u^, 9^, ■•■’®ls^ state variables represent small dis- 

placements and their incremental rates from an equilibrium point. There- 
fore, the minimization of the performance index (2-3) minimizes the energy 
in the state and control displacements. 

The linear optimal quadratic control law yields linear feedback of 
states and results in a reasonable damped, stable closed loop system. The 
quadratic cost function may be regarded as a generalized weighted mean 
square error criterion for the multi -variable system working under the 
penalties on the control variables . 

The quadratic cost fmction for this problem is written with a 
weighting matrix Q for the state (or the output) , and a scalar weighting 
R = for the control as follows : 


J* 


= Min 
u 



[x"^(t) Qx(t) + p^u^ (t)] dt 


(2-4) 
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The optimal control law u(t) is obtained as 


u(t;) = -R"^B'^KxCt) (2-5) 

It has been shmvn that the gain matrix K must satisfy the algebraic 
matrix Riccati equation: 

KA + a"^K - — KBb'^K + Q = 0 ; where Q > 0 (2-6) 

P" 

which yields a feedback control system with the stable, closed loop roots 

R^{X^[F]} H - •^Bb'^K]} < 0 for all i = 1,2, 3, 4 

P 

To solve Eq. (2-6) for K matrix, the weighting matrix Q must be selected. 
For exanple, the longitudinal dynamic equations for a hovering helicopter 
(4th order system), one can write: 

u(t) = -R~^B^Kx(t) = -G^Xj^(t) - G 2 X 2 (t) - G 2 X^(t) - G^x^(t) (2-7) 


where 


K = {k.^} ; i, j = 1,2, 3, 4 

^1 ^^1^21 ^ ^B1^41^ ^ ^ 
^2 ~ ^^B1^22 * ^B1^42^ ^ ^ 
S = <;Sl^23 " ^B1^43^ ‘ P' 
S = C^Bi’^24 ^ ''bi‘^44^ ' P' 
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After the Laplace transform, the control law Eq. (2-7) can be written as: 


Bis(s) = - G^CT^s + 1) x(s) - GgCTgS + 1) e (s) (2-8) 


where 


= '=!■ = "^3> ’’x “ ’'e ' '=4''S (2-9) 

The optimal control system theory attempts to provide an analytical 
design procedure that lessens the designer's load in the design task, and 
locates more of the load on the computational machines. The problem is to find 
an optimal linear controller for a linear plant by minimizing the above 
quadratic system performance criterion or a quadratic cost function. The 
optimal linear-quadratic control theory provides a well organized_design 
procedure for a linear system with a feedback structure. The linear feed- 
back structure would probably be, from the practical point of view, the 
most beneficial result of optimal linear-quadratic control theory. The 
application of the theory, even the multi loop system design, which is often 
difficult in the conventional approach, may be thought to become easily 
accessible for system designers. All the possible feedback loops are taken 
into account and all the corresponding feedback gains are determined in 
the design procedure. The resulting closed loop system appears to be not 
only a stable system, but also a minimum weighted mean square error system. 
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The benefits and the mathematical beauty of the optimal linear- 
quadratic control theory -might give a faise impression that almost all 
the difficulties in linear system design have been eliminated. The 
entire theory has been developed under the presumption that the weighting 
matrices in the quadratic cost function are somehow given a priori. They 
are completely left open for the designer's choice. 

The weighting matrices are usually selected only for their diagonal 
elements by the practical system designer on the basis of his engineering 
experience, coupled with simulation runs for different trial values. One 
can be obliged to do a tediously long trial -and-error approach, especially 
in the design of multiloop systems., Many of the attractive features of 
optimal linear-quadratic control theory seem to disappear.! The optimal 
linear-quadratic control theory may be even worse than the conventional 
methods for multiloop system designs, because the weighting matrices have 
no apparent quantitative relationships with the closed loop system dyna- 
mical characteristics. None the less, the weighting matrices determine, to 
a large extent, on the resulting system dynamical characteristics. Hence, 
the selection of weighting matrices often become a considerably difficult 
task. One may be able to say that most of the potential difficulties of 
fflultiloop system design have been concentrated on the selection of weighting 
matrices, and the optimal control theory does not ease essential difficulties 
in the system design. 

Accordingly, the selection of weighting matrices really has a crucial 
significance in the practical system design through the optimal linear- 
quadratic control theory. Some of the qualitative properties of weighting 
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matrices are summarized in Ref, [26] . They may help designers to select 
the trial weighting matrices . In order for the theory to be practically 
applicable, however, it is necessary to develop a reasonable and unified 
way for the quantitative selection of weighting matrices. This is one of 
the topics considered in this work. 

Two methods for the quantitative selection of these weighting matrices 
have been available, 

1, Method of Authorities of Variables (Refs, \ 2 ^^ and [28]) 

This method is based on the normalization of output variables and 
control variables. Every weighting factor for the output variables and 
the control variables is inversely proportional to the square of the pre- 
scribed maximum for that variable. The quadratic cost function becomes 
a linear combination of normalized mean square errors and controls. This 
method only provides an initial estimate for the weighting matrices, which 
are assumed to be of diagonal forms . However, there are no clear-cut 
reasons for the weighting matrices to be necessarily of a diagonal form, 
since the output components are usually coupled with one another. In order 
to refine the resulting system, a trial-and-error approach would be inevi- 
table. Besides, it is often a difficult task to estimate the maximum values 
of all the output and the control components. Even if it is possible, there 
is no guarantee that all the outputs and the controls of the resulting sys- 
tem are suppressed within the prescribed maximim values. 

2 . Closed Loop Pole Allocation (Refs. [20] and [6]] 

This method achieves the specified closed loop poles by appropriately 
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choosing the weighting matrices in the quadratic cost function. The ori- 
ginal work for this method has probably been done by E.J. Ellert (Ref. [6 ]) 
W.M. Wonham proved in Ref. [35] that it is possible to implement arbitrarily 
assigned poles on the complex plane with some linear feedback control law, 
if the open loop system is completely controllable. However, this has been 
proved independently of quadratic cost functions. More recent studies (Refs 
[35] [37]) on the inverse problem of linear feedback control show that 

every linear feedback control minimizes a quadratic cost function with cer- 
tain weighting matrices. Combining Wonham*s work and the results out of the 
inverse problem, it is possible to place the closed loop poles anywhere on 
the complex plane by choosing appropriate weighting matrices. However, the 
difficulty is how to allocate all (but not some] of the closed loop poles, 
especially in high-order systems. The dominant closed loop pole allocation 
has also been studied (Ref. [20]]. This is a method to approximately obtain 
the desired dominant poles. However, this does not answer the essential 
difficulties in deciding the number of dominant poles, and in choosing the 
dominant poles. Besides, this approach usually yields a system with very 
high feedback gains, and requires complicated computational procedures. 

If the desired closed loop poles can be assigned, the optimal linear-quad- 
ratic control theory may not be necessary, but some other direct approach 
may be more convenient (example: Ref. [31]]. 

Thus, there is still a large gap between the optimal linear-quadratic 
control theory and its application to practical system design. Until such 
a gap is bridged by developing a reasonable, simple, and well-unified method 
to select the weighting matrices, the optimal linear-quadratic control 
theory would remain to be just a "theory" and would not be accessible to the 
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designers of practical systems. 

This work is based upon Ref. [7], which attempts ta bridge the gap. 
Below and in Appendix E is given a summary of the two new methods described 
in Ref. [7] , to select the weighting matrices of general forms 
necessarily diagonal matrices}, 

3 . OWEM Design Method (Ref, [7]} 

Based on the thought that the quadratic cost function is a legitimate 
system performance measure in the linear-quadratic control theory, the 
minimum quadratic cost function with the optimal linear feedback control 
law is actively used as the criterion to select the weighting matrices. 
First, the nature of the minimum quadratic cost function has been examined 
for particular sequences of weighting matrices. It leads to the necessity 
of some constraints on the weighting matrices, as a minimum equal to zero 
occurs when the weighting matrices are zero. In Ref. [7], reasonable con- 
straints are explored to make the optimal weighting matrices not equal to 
zero or finite, and this appears to be a determinant constraint. Under this 
constraint, the minimum quadratic cost function is further minimized with 
respect to weighting matrices. The resulting weighting matrices are based 
upon the system noise characteristics, and hence, the minimum of the quad- 
ratic cost function becomes a function of the noise N or J • CN,Q,R}, The 
resulting closed loop system dynamics is a system that minimizes the effect 
of the noise or disturbance for which it was designed. It was also sho™ 
that the resulting optimal weighting matrices Q and R were inversely pro- 
portional to th.e output covariance matrix, and the induced covariance matrix 
of optimal controls, respectively. Obviously, they are not necessarily dia- 
gonal matrices. Note that these results are quite similar to the weightings 
proposed by Gauss for the weighted least-square estimation problem. 
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The optimal vfeighting matrices appropriately normalize the output 
and the control variables, and also equalize in groupwise all the addi- 
tive terms appearing in the quadratic cost function. The two functions, 
i.e., the normalization and the equalization, are fundamental roles re- 
quired for the weighting matrices in quadratic cost functions. Because 
of the normalizing effect of weighting matrices, the resulting closed 
loop system dynamics became irrelevant to the physical dimensions which 
were used to express the plant dynamical equations. 

It was also shown that the ratio p = |R|/|Q| between the prescribed 
determinants of R and Q represented some kind of measure about the trade- 
off between the contribution of control variables, and that of system out- 
put variables to the quadratic cost function. The tradeoff parameter p 
would play the key role for determining the closed loop system bandwidth. 

The method described above provided a well formulated set of necessary and 
sufficient conditions for reasonable optimal weighting matrices with only 
one unknown parameter or trade-off parameter equal to p. It required the 
a priori knowledge about the system noise intensity matrix N since the 
criterion (N;Q,R) was a function of N. The a priori knowledge about 
the plant noise may not be available in many cases. Then, the second 
method for the selection of the weighting matrices was developed by intro- 
ducing a new criterion which is independent of the noise intensity matrix N. 

The concept of n-dimensional hyper error ellipsoid (state error set] 
is used and a measure of the state error set was defined by its n-diraensional 
volume (V(t)). The "Relative Rate of Convergence of Error Sets" was defined 
as I — ~ lt~ ' ^ ^ represents a sort of generalized system re- 

sponse speed. This new criterion is independent of the plant noise intensity 
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matrix N, and was shoivn to be exactly the same as the so-called "Total 
System Damping". The well formulated necessary and sufficient conditions 
for the noise-free optimal weighting matrices are derived as the solution 
of the Mini -Max problem for the RERACES (or TSD) . 

It has been shown that, using the second method selection of the Q 
matrix yields the same result as the first method, if the noise character- 
istics (see Fig. 1} are assumed to be white noise disturbances with N^=kR 
In other v;ords, a white noise disturbance in u(t). The effect of noise or 
disturbance bandwidths in the helicopter application is discussed in Chap- 
ter 3, and is often small for the range of operating conditions used in 
that application. 

The initial application of this method for determining the Optimal 
Weighting Matrices to second through fifth order example in Ref. [38] in- 
dicates the advantages of this approach. For these examples, the conven- 
tional and optimal methods have been used extensively, and the results are 
well known. For these examples, the area of the pole location for the pre- 
cision hover of a helicopter has been well established, and each of the 
above methods can be applied yielding (often by trial and error), about the 
same results. However, it has been demonstrated, by these examples, that 
the application of the OWEM method yields, without a trial and error method, 
a good and well damped solution to the problem. 

This solution is independent of scaling (dimensions) of the problem 
statement, and assumes in the first instance, that there is no a priori 
different penalty on error, gain, frequency, etc., among the state variables. 
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In siinunary, for the OWEM design method it is required to know: 

1. the system matrix 

2. the control matrix 

3. the noise or disturbance matrix 

The OWEM design method then determines the weighting matrices Q and 
R for minimizing the quadratic performance index under the assun 5 )tion 
that all errors are equally important, and hence, these matrices for the 
calculation of the feedback gains. When the noise N is unknown, Q and R 
can be computed by the mini -maximization of TSD with respect to Q 'and R 
under prescribed determinant constraints. This is equivalent to a white 
noise disturbance with a magnitude inversely proportional to the con- 
trol weighting R, and the minimization of TSD is equivalent to the mini- 
mization of J . (N ; 0,R) . 

min^ e ' 

The OWEM design method is based upon minimization of known or esti- 
mated disturbances and results in a control system with adequate stability 
or damping. In the generalized formulation no consideration is given to 
system limitations, such as acceleration limits, gain limitations, and so 
on. Hence, in the design of a flight controller, the performance and con- 
trol boundaries must be considered. These boundaries are essential in the 
practical application of the design method. For example, take the case of 
a hovering helicopter Clongitudinal) . 

A design is to be made for an automatic control system. However, the 
assumption is made that when the system is switched to manual control, it 
is desirable that the mode frequencies do not change appreciably. Hence, 
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the boundaries for the automatic system should be similar to those of manual 
control. The longitudinal control of the precision hover of a helicopter 
is in principle dominated by 2 modes; a "position" mode corresponding to 
the "outer loop" and a "attitude mode" corresponding to the "inner loop" 
of the feed control system. 

An estimate of the boundaries of the "position" mode is given in 
Ref. [22], and summarized below. 

A massless point controlled according to a linear error control law, 
e = -Ke, would result in a first order system with a closed loop bandwidth 
co^j^ = 1 / t ^^ = K. Instantaneous speed changes not being possible, even the 
simplest model has to incorporate the inertia. This leads to a second order 
system for control. The most significant assunption in this very sinple 
model is that tilting the thrust vector and changing its magnitude, is pos- 
sible with negligible delay compared with the position keeping bandwidth. 

The block diagram of such a simplified model is sho™ in Illustration (2-1 ) . 


Acceleration 



Illustration 2-1. 


Simplified Block Diagram of a Hovering Helicopter 







The input to the longitudinal positioning loop is shown as the 

acceleration disturbance. The' gains gG^ and gG^ determine the accelera- 
tion per foot error and per ft/sec error rate, respectively. The following 
fundamental relationships are easily deduced from the block diagram. 

The natural frequency and therefore, the bandwidth of the position 
loop is 


; C2-10) 

gG is the "spring constant" of the position loop. The relative damping 
ratio is 


?p = Cv^/2) = 2.83 (Cy/G^) C2-113 

The relationship between the bandwidth and the damping ratio is determined 
by 


G*/G^ = 2Cyo)p) (2-12) 

and, for a predetermined bandwidth and damping ratio 
= (1/g) and G- = (2/g) ^ 


In the case of a constant acceleration disturbance, the steady state error 


is ag/w^ = a/G . The sin 5 >le model provides reasonable estimates of the 
P ^ 
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relationship between bandwidth gains and errors as long as the bandwidth 
of the attitude control loops is several times faster than that of the 
station keeping loop . Such a separation is feasible and appears desirable 
in most applications. 

For the choice of bandwidth of the position loop, the following con- 
siderations can be made. A large leads to a high error gain G^, a small 
results in a large steady state error for constant acceleration distur- 
bance. 

The table below shows some characteristic values of bandwidth, gain 
and error. 


Up Crad/secJ 

.1 

.2 

.5 

1.0 

(deg/ft^ 

.018 

.07 

.45 

1.8 

steady state error 





for .Ig acceleration 

320 

80 

13 

3.2 


This table helps to narrow the choice of ( 0 ^ to between approximately 
.3 and .5. With values under .3, the error becomes excessive, above .5 the 
thrust angle variations with even small errors are rapidly becoming too 
large. 

The desirability to make the bandwidth of the "attitude” mode faster 
than that of the "position" mode is limited by the manual attitude control 
loop. This results in an "attitude" mode of frequency in the order of to > 
1.0 to 2 rad/sec. 

In the conventional control system design method, the gain limitations 
are taken into account for choosing the "inner" and "outer" loops. In these 
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hovering helicopter examples a "tight" inner loop (altitude) is designed 
with sufficient damping by adjusting Gg and Gg, the outer loop (position) 
is determined by the position gain limitation G ; the velocity gain G. is 
adjusted to provide sufficient damping. 

In the OWEM design method as developed does not take into considera- 
tion gain limitations. When there are no gain limitations, the operating 
point in the linear quadratic optimal design is determined by the control 
weighting. This number is determined by the control limitations, and the 
maximum exqjected disturbances, and/or command inputs. 

However, the OWEM design method can be easily adapted by additional 
information, into the design. 

Itfhen no additional information is used, the conqjutational procedure 
regards the errors and the associated feedback gains, etc., equally im- 
portant. Hence, in this design method, the root location tends to be cir- 
cular with more damping than the Butterworth configuration. After the 
initial design is made, the resulting mode frequencies, errors, gains, etc., 
can be analyzed. IVhen the value of these quantities are not in conflict 
with the practical design requirements, then one can consider the design 
to be completed. 

However, in many cases, practical physical limitations have to be 
imposed on errors, loop gains, mode frequencies, etc. For example, in 
the hovering helicopter, the position mode frequency is acceleration limited 
e.g., .3 to .5 radians/sec. 

On the other hand, the attitude loop should be a much higher frequency, 
e.g., 1-3 rad/sec. 
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In Appendix E is shown in the calculation of the optimal weighting 

matrices the groups wise equalization matrix (53 is equal to the unity 

T 

matrix. This imposes an equal weighting on the elements of the y* Q y'^. 

In engineering terms»this implies that all the mean square error associated 
control, gams, etc., are of equal importance. 

The initial design is done on the digital computer using S = I. The 
resulting design is then examined with respect to errors, control and gain 
limitations. For example, one can examine errors and mode frequencies as 
a function of the control weighting p. IVhen this examination shows that the 
control weighting p determined by the control limitation, no other boun- 
daries (imposed or desired) are violated, then the design is finished. If, 
however, one finds that for the desired control weighting p, one of the pre- 
determined requirements is not met or gain boundaries are violated, then one 
should determine the violated boundary on the Optimal Root Location (ORL) plot. 

In the example shoivn for the longitudinal hovering helicopter, the 
position acceleration limitation limits the p to approximately equal to 800 
(.3 to .5 rad/sec.), while the desired altitude boundary yields a 
p approximately equal to 3200. This inplies that in order to reach the con- 
trol boundary and acceleration boundary at the same time, one should use an 
equalization matrix which has a value of 1/4 for the state variables associated 
with this limitatioh. 

In other words, one can assign a different weighting on the system 
errors if this also is equivalent to assigning a diagonal Error Scaling Matrix 
(ESM) . The designer is free to investigate the effect of the ESM matrix upon 
the final system design. Note that the ESM matrix is a diagonal matrix, so 
that the elements of the matrix are easily interpreted. 
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In this manner, an objective measure of an error criteria is used 
to determine the weighting matrices . The approach taken in this work 
imposes the boundaries more or less by engineering judgement in terms 
of errors. A very important further theoretical investigation with 
significant practical aspects is the introduction of constraints in the 
OWEM design. 

The OWEM design method is also directly applicable in model following 
for the determination of the weighting matrix Q (Appendix Gl, Eq. Gl-22). 
The weighting matrix Q in the model reference system is composed of the 
original weighting factors Q of linear feedback systems. The OWEM method 
determines the Q matrix for minimization of disturbances. Hence, a weight 
ing matrix determined by this principle should yield a prefilter model 
reference control system that minimizes the errors due to the disturbances 
Another important investigation will be the practical aspects of the 
use of OWEM design criteria for incomplete measurement and feedback of all 
state variables. 
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CHAPTER 3 


CONTROLLER DESIGN FOR A HELICOPTER (OWEM) 

( A Helicopter system and control equations 

In this chapter the optimal control design will be illustrated on the 
controller design of a helicopter. The system to be controlled is a sys- 
tem representative of an uncoupled hovering helicopter CH54B given by Eq. 
3-1. 


r 


•• 








t 

— 

r 

«• - 

* 

u 


X 

-g 

X 


0 

' 0 

0 

0 

r 

0 


u 


u 

q 






) 




« 

9 


0 

0 

1 

1 

I 

t 

0 

' 0 

0 

0 

1 

1 

1 

0 


0 

q 


M 

u 

0 

M 

q 

1 

( 

1 

0 

‘ 0 

0 

0 

1 

t 

1 

0 


q 

* 

w 


0 

0 

0 

1 

\ 

1 

z 

w 

* 0 

0 

0 

i 

i 

0 


w 

- - 



^ » 

« 



- ^ . 

M ^ ^ 

• — « 




- - 

• 

V 


0 

0 

0 

t 

t 

0 

, Y 
' V 

g 

Y 

p 

1 

1 

Y 

r 


V 



0 

0 

0 

> 

1 

0 

' 0 

0 

1 

1 

1 

0 


4> 

« 

p 


0 

0 

0 

1 

1 

1 

0 

' h' 

. V 

0 

L' 

P 

1 

1 

0 


p 


- 

- _ 

. — _ 



"l 

- - - 

■| 

— — 

— 

r 

- - - - 


-- 

r 



0 

0 

1 

0 

1 0 

0 

0 

t 

N-* 

r 


r 


21 




For the optimal design, the following equations are used and solved 
with the aid of a digital computer. 


x(t) = Ax(t) + Bu(t) + Cn(t) 


y(t) = Hx(t) 


C3.3] 


where x(t): n dimensional state vector 

u(t): m dimensional control vector (m<n) 

y(t): p dimensional output vector (p^) 

n(t): i dimensional noise vector (Z^n) 


£{n(t)} = 0 

Cov{n(t ) ;n(x)} = N6 (t-T )( unknown) 
A,B,C and H: nxn, nxm, nx£. and pxn constant 
matrices, respectively. 
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T 


Cost function ; J (NjQ,R) = Min rLim^ f [y^(t)Qjr(t)+u'^(t )Ru(t)]dt} 

uCt) ^< 2 ^ ^ 


-T 


= Tr[CNC^K]. (N : " unkno'wii ) 


If = kC'^BR then 

e 

J ;Q,R) = kTr[BR"^B\] 

THUl 6 


C3.4) 

C3.5) 


where C : Axn pseudo-inverse matrix of C,i.e., 

(=(C^C)‘^C^) 

—1 T 

Optimal control law ; u(t) = - R ^B Kx(t) (3-6) 

where K; uxn unique positive definite solution of ( 5 .7 5 

KA. + - KBR'^B^ + H^QH =0 ( 3 . 7 ) 

The optimal control law (3.7) can be computed without knoA^ing 
H if Q and R are given. 


If the Q and R matrices are estimated, then for these Q and R values the 
optimal control law can be calculated. 

When the values for the matrices Q and R are not estimated, then one needs 
an additional criteria. For example, minimizing total system error in total 
system damping, etc. as discussed in Ref. (7). Assuming that gust or noise 
information is not available, the selection of optimal weighting matrix is 
done by minimizing the criterion L(Q,p^), as discussed in Ref. 7 and Appendix 
E, 


Auxiliary performance index 

Total System Damping = TSD: the negative sum of the closed loop roots 

L-^ = Min Max {TSD} = Min Max {Tr[-(A-BR~^B'^K)] } (3.8) 

Q R Q R 


or equivalently 
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L* = Min Max {Tr[BR~^B'^K] } (’ . * A = constant matrix) 

Q R 

subject to KA + A^K - KBR"^b’^K + h'^QH = 0 
Note; (1) TSD is equivalent to the Relative Rate o£ Change of Error Sets 

CRERACES, see Ref. [7], Chapter IV), and is a generalized system 
response speed. 

[2) The extremum Crainimax) of TSD does not exist unless Q and R are 
appropriately constrained. 


Constraints on Q and R 

iQl = 1, |Ri = 

where 

I CO 1 : "the determinant of (*) 


Optimal weighting matrices 


Q = 



T -1 
CHPH ) 


C3.9) 


(3.10) 


2/m b'^KB - b'^KPKB 
'V'Ib^kb - b^kpkb! 


(3.11) 


where FP + PF^ + BR ^B^ = 0 
F = A - BR“^b'^K 


(3.12) 
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The latter equations are solved on the digital cou^uter for the 
matrix Q. The optimal weighting matrices are then used in the Riccati 
Eq, 3-8 which yields the optimal control law (Eq. 3-7). 

For a hovering helicopter CH54Bj the numerical values for Eq. 3-2 
are shown in Eq. 3-13. 
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These equations can be divided into four separate systems of equati-^ns, 
namely : 

1. Vertical Control System 

2. Yaw Control System 

3. Longitudinal Control System 

4. Lateral Control System 

For a helicopter in the approach phase, the numerical values for 
Eq. 3-2 are shown below in Eq. 3-14. 
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B Controller Design for Vertical and Yaw Control 


In hovering, the helicopters have inherent rate damping in vertical 
and yawing motions with stable vertical speed and yaw rate characteristics. 
However, with regard to the position stability, both the position control equa 
tions are neutrally stable due to pole located at origin. This means that 
the helicopter, if uncontrolled and if exposed to disturbances, will drift 
away from desired hover points. Taking these into account, the following 
design objectives are recommended for precision hovering. 

Cl) To obtain the stable well damped position control system 

C2) Good system response to position command input, with a small or 
no overshoot, and minimizing integral square error (ISE) 

C3) Minimizing root mean square position errors and resulting controls 
to random gust input. 

Some of these requirements oppose each other and a compromise becomes 
necessary. 

The characteristics of the decoupled (uncoupled) motion of vertical and 
yawing equations are of the second order rate command systems as in Eq. (3-14) 

x(t) + ax(t) = b5(t) + (cw (t) ) (3-15) 

& 

where 

x(t) : position z, ft or <j), rad 
x(t) : rate of position 2, ft/sec or r, rad/sec 
a : rate damping Zw or Nr > sec ^ 
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b : control derivatives Zg, ft/sec^rad or Ng , rad/sec^rad 

5(t): control input 5^ or 5^ , rad 

vertical gust or sideway gust, £t/sec 

c : gust derivatives, 1/sec or N^, 1/sec ft 

We assume for calculation purposes that disturbances even in decoupled yaw 
equation has effect on control through N^term. 

Design Method by Conventional Selection of Weighting Factors 

The optimal control is derived by minimizing, in a disturbed system 

the following integral without w 

- S * 

^^11 "" ^22 ^ r6^(t))dt f 3 _ 16 ) 

Where : weighting factor for position 

*^22 * weighting factor for rate 

r ; weighting factor for control 



Selection of Weighting Elements q^^^, q ^2 and r 


In both systems the closed loop system characteristics are covered by 
selection of non-zero diagonal q^^, q ^2 and r only, as shown in Appendix D. 
The weighted quadratic terms in the integral should be non-dimensionalized. 
The following units axe selected for q^^, q 22 


^11 

^22 


Vertical 

Yaw 


ft"^ 

rad"^ 

C3-17) 

ft~^.sec^ 

j-2 2 

rad sec 


rad"^ 

rad'^ 
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The values of weighting factors are specified by choosing an equiva- 
lence of the values of position errors (Z,<j))jrate errors and re- 

sulting control (6^, 6^ ). The selectxon is as follows: 

Vertical [10 ft]^ = (5 ft/ sec) ^ = (0.1 rad)^ (3-18) 

Yawing (0.2 rad)^ = (0.1 rad/sec) ^ = £0.1 rad)^ (3-19) 

These weighting factors correspond to the case (5) with r = 100 in 
Table 1 below, where a number of computed numerical examples are listed. 


Table 5.1. Computed Numerical Examples 
Case No . Weighting factors 



^11 

*^22 

r 

(1) 

O 
1 — 1 

0 

Q Co 

( 2 ) 

1.0 

0.25 

0 ^ 00 

( 3 ) 

o 

I— f 

1.0 

0 -J- OO 

( 4 ) 

1.0 

2.0 

0 -> CO 

( S ) 

1.0 

4.0 

0 00 

( 6 ) 

1.0 

16.0 

0 00 


Note that only the relative magnitudes (values) have effects on the 
characteristics of the closed loop system (see Appendix D) . 


Optimal Control. The optimal control to minimize the performance index (Eq. 
3-16) is given as function of q^^^, and r (Appendix D) . 
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Substituting this into Eq. (3-15) gives the closed loop system 


x(t)+^a^ + 2hyl~^ + ?|ix(t> +b-J"^x(t) = 0 (3-21) 

This derivation and its properties are shown in Appendix D. 

The optimal control, in general, has been solved for the system without 
noise as shown in Appendix D. The separation theory in Refs. ( 20 , 27 > 39) 
warrants the optimal control to be optimal for the system with gaussian white 
noise. However, unfortunately it is not guaranteed to be still optimal for 
the system with nonwhite or nongaussian gust. 

Root Locus for some < 122 '^ *^11 ~ 1*0* root loci with varying r are 

shown in Figs. (6 and 7) for vertical and yaw equations, respectively. The 
complex roots are circvunferently distributed to the origin. These numerical 
examples are of the following numbers . 

Vertical 

a = Z. = -.269 1/sec 

W 

b = Zg = -292 ft/sec^/rad (3-22) 

c 

c = -z = -.269 1/sec 
w 

Yawing 

a = N = -.482 1/sec 
r ' 

b = Ng = -8.95 rad/sec^/rad (3-23) 

°r 

c - -N^ = -.0019 rad/ ft/sec 

Notice that the closed loop system is sufficiently stabilized by weighting 
and r only, and the damping ratio approaches .707 as r becomes smaller. 
These poles are called 2nd order Butterworth poles (Ref. 9 ). Addition of q 22 
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gives rise to a more damped system. The damping ratio and undamped natural 
frequency are given as: 






(3-24) 


C3-25) 


The latter determines the system response speed, or accuracy, which is in- 
versely proportional to quartic root of control weighting r. IVhen the 
weighting q ^2 satisfies the following inequality 


f 


^11 

^22 


Then all roots are on negative real axis, since 


^22 


(3-26) 


has the role 


of a zero in the root locus. It is observed that such selection of weightings 
as in this case 'that is' no cross product terms in Eq. (3-16), place the 
closed loop system poles in certain limited left half plane, A desirable re- 
gion of pole location closely relate with the requirements (2) and (3) . In 
spite of increasing accuracy or speed of response, placing poles too far away 
in the left half plane, by increasing feedback gains or control action, is 
likely to excite the structural vibration mode of the helicopter. On the other 
hand, poles near the origin present a slow response and roughness of control. 
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The feedback gains are shown in Figs . (8) and (9) , which also show 
the corresponding augmented damping and stiffness. The feedback gains 
increase as r becomes smaller, and increase in q ^2 raises the rate feed- 
back gains. 

The optimal feedback controls for both vertical and yawing systems do 
not yield a system which has an overshoot for a step input. Hence, some 
of our design objectives (stable, no overshoot) are definitely obtained. 

ISE and ISU Evaluation for Comnand Input . In general, the ISE and ISU for 
the second order system to a unit command step input are given as 


ISE = 


u) 


c? 


n 



(3- 2 7a) 



x^dt 


x(o) = 1.0 
x(o) = 0.0 


(3-27b) 


ISU 


(0 

n 





(3-28a) 



6^dt 


x(o) = 1.0 
x(o) = 0 


(3-28b) 


where 


^n : closed loop undamped natural frequency rad/sec 
^ : closed loop dan 5 )ing ratio 

a : open loop damping 
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which are sketched as a function of ? in Figure 10 where the minimum ISE 
is attained when S = 0.500, and the ISU is inversely proportional to ? . 

In the s-plane the curves for ISE = const, and ISU == const, are shown in 
Figure (18) . 

Using Eqs. (3-24), (3-25), (3-27) and (3-28), the integrated square 
errors and controls to command step input are computed with varying r, 
some <122*^ “^ll “ 1*0 as shown in Figs. (11) and (12). The ISE value 
approaches the minimum as r goes zero, while the ISU goes to infinity. 

The minimum ISE value of the optimal control system is given by substituting 
Eqs. (3-24) and (3-25) into Eq. (3-27), 


Min . I SE 


1 / ^22 

2 y qil 


(3-29) 


at r = 0 where the ISU becomes infinity. This value would become criterion 
for selection of q^^ and q22' that of the model equation given by 

Eq. (D2-21) in Appendix D. For reference of ISE values, it is desirable to 
compare with the ISE of the first order system of which command response is 
easily depicted in mind. Therefore, the folloiving time constant may become 
an index: 


T = time constant of the first order system with same ISE 

V 

(3-30) 
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which is labeled in Figs. (10), (11), (12), (38) and (63). Similarly , one 
may define the following equation as mean control acceleration. 



which would be a rule of thumb for the upper limit of ISU. It is interesting 
to see that ISE of the optimal control system is a little bit larger than the 
true minimum ISE as shown in Figure 10. 

Gust Response . The requirement for gust response supercedes that for com- 
mand input in position keeping such as hovering. In addition, knowledge of 
the gust response is important since the optimal control is obtained in- 
dependently of the random gust. The assumed gust model is of the following 
power spectrum density {pi.s.d.). 


S 

rr 


( 03 ) 


2g^d 
d^ + 03^’ 


(3-52) 


where 


d : bandwidth of p.s.d. of the random gust 
O' : rms value of random gust 
The numerical examples are given by 
o = 10 ft/sec 
d= .314, 1.0, 2.0 rad/sec 


(3-33) 
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which is a considerable gusty condition. 

The rins position error and control are calculated with varying r as 
shorn in Figs. (13) to (^16) for both vertical and yawing equations. From 
these figures, the gust response would be alleviated by increasing 9 , 22 ^ 
increasing rate feedback (see Figs. 8,9), with less control effort. It is 
much more reduced by decreasing r 'that is' increasing rms control action. 

As r goes to zero, the rms control force (acceleration) increases to cancel 
the random gust force with little excess as shown in Figs. (13) and (14). 
Note that excess control action is caused by control system time lag. 

For different bandwidths of p.s.d. of random gust under its constant 
intensity, the rms position errors and controls are studied in Figs. (15) 
and Ci6). These results show that random gust with higher frequency compo- 
nent bring less effect on both rms position error and resulting control. 

These results are explained by the frequency domain analysis as follows: 
The rms value of the gust response is described as (Appendix Gl-3) 


RMS value = 




iGCjoj) 


S^^(w) du 


(3-34) 


where G{j(o) is the frequency response of the transfer function to the gust. 
For example, the frequency responses of the transfer function of ver- 
tical position error and control to the gust input are shown in Figs. (17) 
and (18). The following three cases are pictured here: 
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Weighting factors 


Undamped 

Case No. 

^11 

^22 

r 

Damping 

Natural Frequency 

CD 

1.0 

0 

13130 

0.71 

0.9 

C2) 

1.0 

0 

8190 

0.71 

1.8 (3-35) 

(3) 

1.0 

4.0 

8190 

1.9 

1.8 

The broken 

lines 

indicate 

the p .s .d. of the 

random gust. The rms position 


error in case (1) becomes larger than the others because of the higher gain 
as shown in Fig. (17) . On the other hand, the control action becomes smaller 
since the relative magnitude of control is smaller. In case (3), though the 
control response can be influenced by a high frequency component of random 
gust, the gust under consideration does not include such frequency component, 
and therefore the rms resulting control is smaller compared to case (2). 

Similarly, it is recognized that the random gust with higher bandwidth of its 
power spectrum density has less effect on both rms position errors and con- 
trol. 

Now, assume that requirements are chosen as follows : 

rms vertical position error ^ 5 ft (3-36a) 

Integral squared position error < 5 ft^ sec(T =1.0) C3-36b) 

— 0 

Integral squared control ;^9.0 ft^ sec~^ (Uni=3.0 ft/sec^) (3-36c) 
Eq. (3-36a) would supercede the others in keeping position. 

The lower limit of r is prescribed by only ISE requirement. We should 
choose q^ij^ and q 22 necessary, 

q22/^llll-0 (3-37) 
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and r should satisfy the following inequality 


1200(8.5)^ — <; 2,200(6.2} (3-38} 

^11 

for = 1.0 and q ^2 “ 0.0. 

The requirements for yawing cases are chosen as follotirs : 

rms yawing error < degrees (3-39a) 

ISE <.50 deg^sec (T =1.0 sec) (3-39b) 

ISU £9. deg^sec"^ (%=30 deg/sec ) (3-39c) 

Then, we should choose as necessary 


qj/lllH-O C3-40) 

and r should satisfy, for q^^^^ = 1.0, q 2 ^ = 0 

1.1 (8.5) ir/q^l 13(5.2) (3-41) 

Note that numerical numbers in the parenthesis indicate ratios of expected 
controls and position errors due to weightings, and the ratios are almost the 
same in both cases. There is a big difference between selected r due to the 
differences in control derivatives. 

The desirable root location to meet all requirements (Eq, 3-36), which 
varies as requirements alternate, is sketched in Fig. (18). 

Owem Design Method . The optimal control is given by, as shown in Appendix 
E2. 



(t) 




5 

1 2 
b3 t3 


) i(t) 


(3-42) 
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Note that the optimal control is determined only by r which is called 'trade 
o££ parameter' in the Appendix E. This optimal control is also designed to 
be optimal for the system with a gaussian white noise as sho\m in the Appendix 
E, but not for a nonwhite gust. 

The close loop system is given by, (See Eq,E2- ) 


x(t) + 



4 * 


3b 


4 

3 


2 

^3 


x(t) 


r -2 
3 


x(t) = 0 


(3-433 


of which undamped natural frequency and damping ratio are: 


w = 
n 



(3-44) 



Root Locus for Varying r 

The root loci are's'ketched in broken lines for vertical and yawing 
equation in Figs. (6) and (7), respectively. The closed loop system app- 
roaches the system with damping ratio = .866 while the poles are pushed away 
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in the left half s plane as r goes to zero. The dxstance of the complex 
roots from the origin is inversely proportional to cubic root of control 
vveighting r. By this method, no overshoot and sufficiently stable system 
can be obtained and some of our design objectives are achieved. 

ISE and ISU for Command Input 

The ISE and ISU values are shoivn in broken lines in Figs. (11) and (12), 
respectively. These(lSE and ISU with varying r)are calculated by substituting 
(3-44) and (3-45) into (3-27) and (3-28) and as r goes zero, obviously 

ISE „ = 0 (3-46) 

r=0 

= infinity (3-47) 

The requirements described in Eqs . (3-36b,c) and (3-39 b,c) would be satisfied 
by choosing r properly. 

Gust Response 

The rms position errors and controls of the gust responses are also 
shown with broken lines in Figs. (13) and (14). Decreasing a weight r on 
control or increase in available control reduces the rms position error. 

The requirements described in Eqs. (3-36) and (3-39) are satisfied 
by choosing 


2700(5.6)5r£7500(3.4) 

2.6(5.6)<r<5.8(3.5) 


(3-48) 

(3-49) 



Note that the numbers in parenthesis indicate the ratio of expected controls 
and state errors due to r called tradeoff parameter in the single input cases 
[see Appendix E] . 

In this way, all the design requirements would be satisfied by this 
method. 

Helicopter in Approach Phase 

The numerical values for the helicopter in the approach phase (Eq. 3-14), 
are not sufficiently different from those in Eq. 3-22 and 3-23 to employ a 
control system with variable gain, especially as the control system is of the 


first order. 



C Controller Design for Longitudinal Equation 


A hovering helicopter can have an unstable, oscillatory mode in 
the longitudinal equation, as shown in Figure 2. Furthermore, the position 
control equation has another pole at origin. That places much workload on the 
pilot for precision hovering, especially in gusty conditions. Therefore, the 
design objectives become, for the automatic controller: 

1. To obtain the stable system 

2. Good response, such as little-overshoot, minimizing the ISE, IS0, 

ISU to a command input 

3. To minimize RMS of the position error, attitude error, and resulting 
control to gust input 

The longitudinal equation is given by 


u_xu_x 0 + g 0 = X„, B, -XU 
u q * Bis Is u g 


- M u+ 0 _ M e = M„, B, - MU 
u q Bis Is u g 


In the state variable form, including position variables 


C3-50) 


r* * n 
X 


0 

1 

0 

0 


X 


~ 0 " 


'o 

• 

u 


0 

X 

u 

-g 

X 

q 


u 

+ 

^Bls 

Is 

-X 

u 

• 

0 


0 

0 

0 

1 


0 


0 


0 



0 

M 

u 

0 

M 




M 

Bis 


-M 

u 


u 


o- 

a 


(3-51) 
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where 


Is 


X 

u 

0 

q 

B 

Xu 

Xq 

Mu 

Mq 

X 


Bis 


M. 


Bis 


longitudinal position 
longitudinal velocity 
pitch attitude angle 
pitch rate 

longitudinal cyclic pitch 
gust input 


ft 

ft/ sec 
rad 

rad/ sec 
rad 

ft/sec 


rate of change of longitudinal force, 1/sec 
with velocity 

rate of change of longitudinal force, ft/sec 
with pitch rate 


rate of change of pitching moment 
with velocity 

pitch rate damping 

longitudinal control derivatives 

longitudinal control derivatives 

gravitational force 


, rad/£t sec 
, 1/sec 

2 

, ft/sec /rad 
2 

, rad/ sec /rad 
, 32.2 ft/sec^ 


The performance index to be minimized is given as follows; 


CqiiX^ q22U^ + q44§^ + dt (3-52) 

Selection of q^j|_Ci = 1>2,3,4] and r. 

The selection of diagonal Q or q^^ places the poles of the closed loop 
system in a restricted part of the left half s plane. (Appendix E). The 
following tinits are selected for q^^^ , q 33 > q 44 ^d r. 
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'^11 

ft~^ 

*^22 

ft'^sec^ 

-2 

*^33 

rad 

,-2 

^44 

rad sec 

r 

rad~^ 


Their values are specified hy choosing an equivalence of x, u, 0, q and 
in the integral of Eq. (3-52). The selection is as follows: 

1 ft = 1 ft/sec = 1 degree = 1 degree/sec =0.1 rad 

(3-54) 

Equation (3-54) is expressed as Case No, (7) with r = 100 in Table 3.2 where 
the selected diagonal weighting q^^^ (i = 1,2, 3, 4) are shown. Note that attitude 
errors 6, 6 are specified in degrees , while in the equation the units are radians. 


Table 3-2. Numerical Examples for q. . 


Case No. 

^11 

^22 

^33 

^44 

(13 

1.0 

— 


— 

(23 

1.0 

— 

820 

— 

(33 

1.0 

— 

3,280 

— 

(43 

0.1 

— 

3,280 

— 

(53 

0.01 

— 

3,280 

— 

(63 

1,0 

— 

3,280 

3,280 

(73 

1.0 

1.0 

3,280 

3,280 

(8) 

0.1 

0.1 

3,280 

— 

(9) 

0.1 

0.1 

3,280 

3,280 

(10) 

0.1 

0.4 

3,280 

13,130 
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Root Square Locus. 


The root loci for the longitudinal system is shown in Figs. Cl 9) through 
{21') . A stable system is obtained using only a weighting of the Q matrix 
and r. Two of the poles approach the asymptote C = 0.707 while attitude poles 
approach the zeroes, as shown in Fig. [19). The root loci starts from stable 
poles which are then shifted to mirror image in the left half plane of the origi- 
nal unstable poles. Increasing the attitude weighting factors from a zero 
value, induces the more damping to the attitude poles, and further increases in 
q ^2 moves the attitude poles away from the asymptote and shifts the position 
poles into the left half plane as shoim in Figs. (20) to (23). As a result, the 
attitude control or inner loop control becomes much faster than the outer loop 
or position control. The boundary condition, whether or not the locus of the 
attitude poles bends inwards, depends on the weighting ratio of q^^^^ and q^^- 
The zero locus is dra^m as a function of in Fig. (22). This indicates 

that an increase in the relative weighting q^^ to q^^^ brings two zeroes to the 
origin* Addition of either or q^^ yields another zero on real axis and pro- 
vides more damping to both the attitude and position poles, as shown in Figs. 

(24) to (27) . Especially, imposition of q^^ 0 locate a zero near 

Hence, only one root goes to infinity along the negative real axis, while the 
others approach the zeroes, and the system is more damped. The closed loop sys- 
tem with damping ratio more than 0.707 would be achieved by the adequate selec- 
tion of q 22 » *^ 33 » *^44 where obviously more weighting would be imposed 

on <^33 

By the OlVEM method, the root locus is uniquely determined with choice of r 
only. Four poles move to infinity along ? ss .65 and C. « .95, and the closed loop 
systems are sufficiently damped. 
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In this some of the design requirements are attained without diffi- 

culty. The feedback gains are shoivn in Figs. (28) to C^7) . Using these 
diagrams, the optimal control is obtained as: 


*1S ■ V * V * * V 


(3-55) 


At infinite r, both feedback gains Kg and approaches nonzero gains, 
so that the unstable poles are shifted to the mirror image in the left half 
plane. For the conventional selection method of and r, the outer loop, 
feedback gain is given by: 


K 

X 


^11 


(3-56) 


which specifies the systems accuracy. From this point of view, Eq. (3-56) 
helps select q^^ and r. 


ISE, IS9 and ISU for a Command Input . 

The integral squared position and attitude error to a 10 ft. command 
input is shoi-m in Fig, (38). In position control, the hovering helicopter is 
desirable to follow the command input as fast as possible without excessive 
attitude change or its rate. As r goes to zero, the ISE value approaches the 
minimm given by that of model equation corresponding the zeroes, while the 
ISO reaches the maximum as shoim in Fig. (38). These values at null r are 
sketched as a function of ratio of weighting factors q^^^ and q^^ i^i Fig. (39). 
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They show that the reduction in ISE increases IS6 in Fig. C29)* Hence, a 
compromise is necessary. The addition o£ the weighting factor, de- 
creases the IS6 considerably, with little increase in the ISE. Finally, 
the amount of control available in Fig. (40) is limited by the power limita- 
tion and mechanical strength of the helicopter. 

As shown in Fig. (38), the integral squared control increases as r 
decreases and decreases more when a weighting is placed on the attitude or 
its rate. 

It is observed that the OWEM has features similar to, but better than. 
Case 1 in Table 2, indicating ISE = 0, ISO = ISU = <n at r = 0. 


Gust Response 


The rms position and attitude errors to the random gust are shown in 
Figs. (41) and (42), for some and varying r. The assumed gust model is 

of the same as Eq. (3-33) and numerical example is made assuming quite strong 
gusty conditions. 


a = 20 ft/sec 
d = .314, 1.0, 2.0 rad/sec 


(3-58) 


The rms position error decreases as t/q^^ becomes smaller and depends slightly 
on Effects of and q^^ are negligible, and hence, not shown in the 

figure. On the other hand, the rms attitude error, due to the gust, depends 
strongly on ^.nd selected q^^^^ (i = 1,2, 3, 4). Adding a weighting factor, q 22 

or q^^ improves the error considerably. In the range of q^^ considered here, 
the ratio of <l33/^jl^2 entirely the order of magnitude of the rms value. 
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The rms longitudinal cyclic control and position error is shorn in 
Fig, (41^. The root mean squared position, attitude errors, and different 
bandwidths of the power spectrum density for a gust with constant intensity 
are computed. It is clear that the rms position errors are less influenced 
by the gust with higher bandwidth. On the other hand, the rms attitude errors 
are affected less at a large r, but become more affected at smaller r by the 
gust components of higher frequency, as shown in Fig. (67). 

The gust responses for the OWEM, shown with broken lines in Figs. (41) 
and (42) are better than Case 1. 

A ranking of all requirements for precision hover can be done as shown 
below. A compromise between them must be found as it is desirable that all 
requirements are simultaneously satisfied. Take, for example, the following 
case: 


root mean squared position error 
root mean squared attitude error 
integral squared position error 
integral squared attitude 
integral squared control 

acceleration limits 



1.0 ft. 

(3-57a) 

< 

0.3 deg. 

(3-57b) 

< 

400 ft^sec (xe<8) 

(3-57c) 

< 

16 deg^sec (6m<2 deg.) 

(3-57d) 

< 

0.05 rad^sec rad) 




(3-57e) 

< 

none 



and (3 -5 7b) vary depending on rms 


The right hands ide of Eqs. (3 -5 7a) 
gust velocity or its rate. The upper limits on Eqs. (3-57c) to (3-S7e) are 
altered proportionally to squared magnitudes of command inputs. However, 
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T , 0 , and U may be independently prescribed, 

max 

Equation (3-57a) exclusively determines the upper limits of r, though 
ISE and rms attitude error requirements may influence this considerably. The 
lower limit of r is described by Eqs. (3-57d) and C3-!i57e-}-. The way of choosing 
*^ii found out by requirements (3-57b), (3-57d) , and (3-57e). 

If one considers the above requirements, then the following set of q^^ 
and r can be chosen as one of desirable sets of q^^ and r. 

r = 2.0 • lO'" 


qil .1, qg2 - 3280 


(3-58) 


With the feedback gains 
K =-.0071 

X 

= -.026 (3-59) 

Kq = 1.58 
K = .56 

q 

It is easily understood in the process of choosing these sets that Case 
(5) does not satisfy Eq. (3-57c) and Case (1) is unsatisfactory_with Eq. (3-57d), 
because of excessive attitude error. An OWEM method which does not take into 
account the difference in position and attitude error weighting should also 
be unsatisfactory (for example, see Eq. (3-S7d) . 

Helicopter in Approach Phase 

The numerical values for the helicopter in the approach phase are given 
in Eq. (3-14). The root location for the optimal roots as a function of the 
control weighting r are shown in Figure 67. It appears that the difference 
in numerical values compared to the hovering case is not sufficiently large 
enough to justify a control system with variable gains. 
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OWEM Design Method 


Applying the OlVEM design method for the longitudinal motion of the 

helicopter gives a root location as a function of r = p obtained by a 

computer solution as sho™ in Figure 19. 

In the initial approach of an OWEM design, all the constants of 

T 

the group-wise equalization of the terms of {y* Q y*} are chosen to be 
equal. This corresponds to an equalization matrix S = I. Making the 
equalization matrix not equal to the identity matrix corresponds to an 
equalization factor of the diagonal terms of the matrix q^^ . This prin- 
ciple is illustrated in this example in Figure 19. 

In this figure is shown the root location obtained by the initial 
application of the OWEM method (equalization matrix = identity matrix). 

The desired root location with associated optimal Q matrix is determined 
by the acceleration limitation of the helicopter (corresponding to an tu 
of .3 to 5 rad/sec). Take, for exan 5 )le, that an equalization was desired 

of one to two for the ratio of the position error to attitude error. The 

• • 

diagonal terms of the Q matrix are associated with the variables x, x, 0, 0 
respectively. Hence, the equalization is chosen to be 1:1:4:1. 

Figure 69 is shown, the root location obtained by the OWEM method, modi- 
fied with an error scaling matrix (ESM) of 1:1:4:1. For this ESM Q matrix a 
root square locus diagram is given as a function of the control weighting p. 
The ESM matrix should also change the weighting of the derivatives of 9, as 
this increases the damping. By trial and error, one can investigate the error 
scaling matrix for the 9 error weighting. In this particular case, 1:1: 4:4 
was chosen. 


49 



In the same illustration is shown the chosen optimal Q matrix, modified 
with an equalization matrix of 1:1: 4:4. 

For this equa-Hzed Q matrix a root square locus diagram is given as a 
function of the control weighting r = p. (curve no. 2) 

For a control weighting of 3200, the (negative) feedback gains are: 

= -.0015 rad/ ft = -.084 deg/ ft. 

= -.0074 rad/£t = -.420 deg/ft. 

Gg = .56 rad/rad = .56 deg/deg. (3-60) 

G^ = .42 rad/rad = .42 deg/ deg. 

The root locations are shown in Figure 69, with the associated feed- 
back gains in Figure 72. The corresponding error and control measures 
are shown in Figures 40, 41 and 42. 

It should be emphasized that, in this work the acceleration limitations 
have been introduced as an "afterthought” in the OWEM design. In the OWEM 
calculations the acceleration limitations were simulated by an "error 
weighting". In principle, it should be possible to introduce constraints di- 
rectly into the OlVEM design method, thereby eliminating the trial and error 
aspect. 

In the hovering helicopter example, feedback gain limitations on the po- 
sition state variable (equivalent to an acceleration limit) should more or less 
"freeze" the position mode when a certain control power is reached. A fur- 
ther increase in control power should only affect the 6 mode. This procedure 
would eliminate the trial and error procedure of the 9 error weighting. How- 
ever, these theoretical aspects need first to be further investigated. 
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D Controller Design for Lateral Equation 


In a hovering helicopter, the lateral equation uncoupled with yaw be- 
comes unstable and poses a comparatively short period mode (see Figs. (2,33). 
Therefore, our design objectives are composed of the same ones as the longi- 
tudinal case. The lateral equations are given by; 


v-Yv-Yp-g<ii = Y., A, -Yv 
V p-^ Als Is V g 


- L V + $- L(p = L., A, - Lv 
V ^ p^ Als Is V g 


( 3 - 61 ) 


In the state variable form, including position variable 
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where 


y : lateral position 

V : lateral velocity 

({) : roll angle 

p : roll rate 


,ft 

, ft/sec 

,rad 

,radsec 
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Als 


Als 


g 


lateral cyclic pitch 
gust input 

rate of change of lateral force 
with velocity 

rate of change of lateral moment 

with roll velocity 

rate of change of lateral moment 

with velocity 

roll rate damping 

lateral control derivatives, 

lateral control derivatives, 

gravitational force 


,rad 
, ft/sec 

-1 

,sec 

-1 

,sec 

,rad/sec ft 
-1 

,sec 

_2 

, ft. sec rad 

T / 2 

,1/sec 
,32.2 ftsec 


The performance index to be minimized is given by: 

+ q22^^ 

Selection of i = 1,2, 3, 4 and r. 

The diagonal weighting matrices Q might place the poles of the closed 
loop systems in restricted parts of the left half plane. The selection of 
Q would be done in the same way as in the longitudinal case and the same 
q. .'s are used as numerical examples [see Table 3.2]. 
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Root Square Locus 


The root loci are shown in Figs . (43) to (51) . The stable closed loop 
system is obtained by a weighting factor only. Addition of gives 
more damping to the system, especially to attitude control loop. Further 
increase in finally moves the poles of the attitude loop further in the 
left half plane, while the poles of position move to the right, 'that is’ 
the position loop reduces stability as shown in Figs. (45 to (47). It is 
seen in Fig. (46) that the attitude poles start going far away for a set of 
q^j^=.l, q22=3280 unlike the longitudinal case. These tendancies seem de- 
sirable since the inner loop control is desirable to be much faster than the 
outer loop position control. The broken line of Fig. (44) indicates the 
zero locus as a function of the parameter of 153/*^ 3^2' starts from the 

point ?=-2.33±j4.06 and approaches the zero along the asymptote ?=.707. Addi- 
tion of (\^2 yields another zero on negative real axis and provides more 

damping to both the attitude and position poles as shorn in Figs. (48, 49, 50 
and 51). Selection of q.,- locates the zero in the far left of the negative real 

Pi? 

axis. On the other hand, selection of q.. locates the zero on =-/ . As a 

^44 1 ^q^4 

result, only one root goes to infinity on the negative real axis while the 
others approach zeroes and the system is more dan^jed. The closed loop system 
with damping ratio more than 0.707, if desired, would be achieved by adequate 
selection of <\ 22 ’ *^33’ *^44 


The feedback gains are shown in Figs. (52) to (61), where it is noted 
that the augmented dancings are equivalent to those of longitudinal ones at 
the same though feedback gains being smaller than longitudinal ones 
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except Ky. given by 


K 


y 



C3-64) 


which specify the system control precision, 
nonzero simularly as in longitudinal cases, 
the optimal control is obtained as : 


At infinite r, K. and K are 
<P P 

Using these gain diagrams , 



(3-65) 


ISE, IS4>, and ISU for Command Input 

The integral squared position error, attitude and control to a 10 ft. 
command input is shown in Figs. (62) to (63). Effects of r and q^^ on these 
results are almost the same as the longitudinal case. When compared to the 
longitudinal case, the IS9 came out larger in spite of the use of less con- 
trol amount (ISU) due to the larger control derivative. Accordingly, the 
command input response becomes faster, that is, less ISE. The min. ISE and 
max. IS0 at zero r are nearly equivalent to those of the longitudinal case 
because both model equations are almost the same. Note that the ISE goes 
zero while the IS0 increases to infinity as r goes zero for the OWEM. 

Gust Responses 

Figs. (64) to (67) show that rms lateral position error (ISE), attitude 
error (IStf*) and resulting control (ISU). From these figures, it is confirmed 
that the lateral motion is much stronger influenced by gust through large 
gust derivatives, and L^, compared to longitudinal motion. The control 
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necessary to trim the translational force and upset moment due to gust 
amounts is twice as large as that of a longitudinal case. Keeping lateral 
posxtion as precise as longitudinal control requires more feedback gains 
corresponding to smaller r. 

Figure (66) indicates the effect of gust with different frequency 
components under its constant intensity to attitude errors. It turns out 
that the gust with cut-off frequency d = 1,0 rad/s influences the lateral 
attitude most, while the pitch attitude is most influenced with cut-off 
frequency d = .314 rad/s at large r. These characteristics are explained by 
the fact that the basic helicopter dynamics have modes excited around fre- 
quency of .36 rad, in longitudinal motion and frequency of .85 rad. in la- 
teral axis. 

Take, for example, the following case similarly to the Eqs. (3-57a) to 
(3-57e) : 


rms position error 
rms roll attitude error 


rSE 

IS9 

ISU 


Acceleration limits 


< 

< 

< 


1.0 ft. 


<_ 0.3 deg. 


(3-6 6a) 
(3-6 6b) 


400 ft sec (T <8 sec) 

2 (3-6 6c) 

16 deg sec (0 £ 2 deg) 

(3-66d) 

0.05 rad sec (.\. < .05 rad) 

~ (3-6 6e) 

none 


The design approach is similar to that of the longitudinal case. 

For the closed loop system, to satisfy all requirements, the following set of 
q. . and r is chosen as one of the desirable sets of q. and r 
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C3-67) 


r = 1.5 • 10^ 

qil = .1, q33 = 3280 

Then, the feedback gains are: 


K - 026 rad/ ft 

X 

=-.087 rad/ ft/ s 
K, =-4.8 rad/rad 

4 > 

K. =-.46 rad/ rad/ s 

<P 


(3-68) 


In a practical application, the acceleration limits will show that 
these gains are too high. 

Helicopter in Approach Phase 

The numerical values for the helicopter in the approach phase are given 


in Eq. (3-14) . The root location for the optimal roots as a function of the 
control weighting p is shown in Figure 68. It appears that the difference in 
numerical values, compared to the hovering case, is not sufficiently large 
enough to justify a control system with variable gains. 

QWEM Design Method 

Using the approach similar to that of the longitudinal case, the root 
locations using the OWEM design is shown in Figure 43 . For an error scaling 
matrix, chosen was 1:1:4:3, and, the feedback gains for p = 1300 are: 

K =“.00075 

y 


(3-69) 


K =-.0025 

y 

K. =“.18 

<P 

Kr =“.09 

<P 

The root locations as a function of r = p are shown in Figure 70 with the 
associated feedback gains in Figure 72. The corresponding errors and con- 


56 



trol measures are sho\m in Figures 62, 63 and 64. 


57 



CHAPTER 4 


CONCLUSIONS AND RECOMMENDATIONS 
The optimal linear-quadratic control theory provides a well 
organized design procedure for the automatic controller design of a 
multivariable system. In this design process an optimal linear 
feedback control law is found, which minimizes a quadratic cost function. 

It is usually assumed that the weighting matrices in quadratic cost 
functions are given a priori. However, the resulting optimal feedback 
control system is a function of the selected weighting matrices. The 
functional relationships between the weighting matrices and dynamical 
characteristics of the resulting system are rather obscure. Therefore, 
the selection of weighting matrices is often a major problem in the 
practical applications of the optimal linear-quadratic control theory. 

In this report the weighting matrices are selected by the OWEM 
design method as described in Ref. [7] and Appendix E. Also in this work 
the weighting matrices using the OWEM design method are obtained by direct 
calculation and not mainly by a creative trial and error approach. The 
OWEM design method is con^ared with other methods for finding the weighting 
matrices . The results indicate that the use of this approach gives at 
least as good results as obtained by other classical design methods. 

Much work still has to be done, especially in the case of multivariable 
input-output systems, and a priori imposed limitations of the practical system. 
However, from the examples in this work, it can be recommended that the OWEM 
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design method is a good candidate to use for the selection of weighting 
matrices. Moreover, the OWEM design method does not require any lengthy 
and laborous estimations of weighting matrices in quadratic cost function. 

A completely computerized system design becomes possible, which can provide 
quick results even for complicated systems. The decision which has to be 
made by the designer after the computerized system design is to choose an 
appropriate trade-off parameter by considering the corresponding feedback 
gains, system dynamical characteristics and other practical situations. 
Practically, of course, further refinements for such a computerized system 
design would have to be done for a few choices of trade-off parameters. In 
the long run, when the given system is highly complicated, for example, a 
system with many feedback loops and/or multiple controls, the computerized 
system design procedure associated with the OWEM method can save large amounts 
of time usually required in the early stage of system designs . 

In conclusion, the OWEM method provides an objective measure to determine 
the weighting factors in the performance index based upon an error criteria. 
Initially, the system errors are assumed to be equally important, however, di- 
rectly from missions requirements one can determine which are the important 
errors, and use this additional information together with the other system re- 
quirements to obtain an improved design. The resulting OWEM design is a con- 
trol system that minimizes the errors due to disturbances and has good transient 
responses . 

This work explores the application of the OlVEM design method to an auto- 
matic controller design for a helicopter at hover and approach. The examples 
demonstrate that this approach is feasible and is a good candidate for the 
control system design. 
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Areas of further studies are: 


1. Further investigation of the practical aspects of system dis- 
turbances in the selection of the weighting matrices, i.e., how 
are the weighting matrices and the resulting optimal feedback 
gains changed if the noise characteristics are varied. 

2. Investigation of the effect of a priori imposed limitations 
[control and/or system boundaries) on the OWEM design method, i.e., 
in this work imposed limitations were combined with an assigned 
error scaling matrix ESM derived from physical limitations. Fur- 
ther investigation of the practical aspects of a priori imposing 
design limitations would be highly desirable. 

3. An investigation of the practical aspects of the effects of an 
incomplete measurement system and feedback of the state variables 
in the OWEM design method. 

4. Further investigation of the application of the OWEM design method 
in model following control systems. In theory, the OWEM method can 
yield the weighting or Q matrix for the prefilter model reference 
control system. However, the practical aspects should be investi- 
gated. 
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Figure 1. Block Diagram for Decoupling 

System dynamics; x(t} = AxCt) + Buft) + Cn(t) 

yCt) = Hx(t) 

where x(t) : n dimensional state vector 

u(t) ; m dimensional control vector 
yit) : p dimensional output vector Cp^n) 
nCt) ; H dimensional noise vector 
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APPENDIX Ai 


review of helicopter dynamics 

Introduction 

In a synthesis or analysis of a control system of helicopters, we are 
confronted with the complexity of the overall helicopter dynamics. There- 
fore, sin5)lifications often are made and are considered the simplified equa- 
tions. However, this can be an oversimplification and when the system is 
constructed, the resulting system may have undesirable properties. 

Although the higher order equation can be manipulated by the computer, 
the analysis and synthesis cannot be easily done. Also the physical in- 
sight into the dynamics is often lost during a process of computation. Fin- 
ally, we would be obliged to use brute forces of computing all parameter 
variations in order to know their effects. The synthesized controller may 
become complex, and often becomes more expensive and less reliable. In prac- 
tice, the precision of hovering can be diminished because of coupling ef- 
fects . 

In general, the complicated dynamics are due to many coupling terms 
which come from asymmetrical body and power sources. With the exception of 
a few terms such as fL and N ) which contribute to dynamic stability, they 

'■XV 

usually bring undesirable effects into the dynamics. These effects can be 
shown from a sensitivity analysis on coupling terms. 

In order to circumvent these things, the dynamics are required to be 
decoupled. The control of the helicopter then becomes simpler and easier. 

The analysis and synthesis can be applied to such conventional techniques 
such as root locus and frequency methods which basically deal with the single 
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input-single output system. The dynamics can be evaluated by hand calcula- 
tion. Even in such modem control techniques as optimal control using state 
space concepts , the analysis and synthesis technique becomes more tractable . 
This method will be dealt with in the other appendices. 

Helicopter Dynamics 

In general^ the helicopter dynamics are \^ritten with respect to the sta- 
bility axis [Ref. Al] . The following small angle conditions are assumed and 
higher order terras are neglected. 

sin 0 = 0, cos 9 =1, sin <}> = ']’> cos cj> =1 , (Al-1) 

0 = q>^ = p, t = r 

u + gcosY^e = X u + X V + X w + X p + X q + X r + AX 
0 u V w p-^ q^ r c 

V - gcosY^fj) = Yu+Yv+Yw + Yp + Yq + fY-U)r + AY 
o^ u V w p-^^ q^ ^ r o*^ c 

w + gsinY 9 = Zu+Zv+Zw+Zp-i- fZ + u;)q + Zr+AZ 
0 u V w p^ ^ q o ^ r c 

. I 

• * xz 

p - r -= — = LutLv+Lw+Lp+Lq + Lr + AL 
u V w p'^ q^ r c 

q = Mu+Mv+Mw+Mp+Mq+M +AM 

u V w p-^ q^ r c 

. I 

r - p ■= — = Nu + Nv + Nw + Np + Nq+N +AN 
^ u V w p^ q^ r c 

CAl-2) 
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where 


AX„ = B, + X. A, + X. 6 + X^ 6 
c I’s A- Is 6 c 6 r 

is Is c r 


AY = B, + Y, A, + Y. 6 + Y- 6 
c B, Is A, Is 0 c r 

Xs Is C 


AZ = Z„ B, + Z. A, + Z. S+ Zs 6 
c B, Is A, is o_ c 


is 


Is 


r r 


AL = L_ B, + L. A^ i)s; h'r— 5 

c B, Is A, Is c I 

Is Is c 


AM = B, + M, A, + Mr 5^+ Nfe 6 
c B, Is A, Is c «r i 


Is 


Is 


AN = N_ B, + N, A, + Nr 6+1^ 6 

c B- Is A, Is 0 c r r 
Is Is c * 


The new derivatives for cancelling the inertia coupling terms 
duced into the last two equations of (Al-2)[Ref. AZ] . 


-p + L'^u + L"v + L^w + L^’p + L'^q + L^r + AL'' = 0 

^ u V w p^ q^ r c 

-f + N'^u + N'*v + N^w + N'^p + N'^q + N'^r + AN^ = 0 

u V w p^ q^ r c 


where 


L = 
u 


L .^N 

U Iv U 


1 - 


CL 

^ XZ'^ 

I I 
X z 


L = 


L + N 
^ V 


1 - 


«xP^ 

I I 
X z 


• • « 
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CAl-3) 


are intro- 


(Al-4) 


(Al-5] 



N" = 
u 


N + ^ L 
u I7 u 

a~i^ 
, ^ xz 


I I 

X z 


N 


xz 


N" 

V 


V 

1 


Jz V 

CX 

^ xz 

I I 
X z 


. . etc , 


In the same way. 


AL" = L" B, + l ; a, + I r 6 + l; 0 
c B, Is A-,^ Is o c or 

IS JlS c T 


AN" = N" + n: B, + N 6 + N; 6 

c B, Is A Is 5 c 5 r 
is is c r 


(Al-63 


where 


xz 


N. 


Is 


Is 


Is 


1 - 


CW 

I I 
X z 


0 




I I 
X z 


.. etc. 


(Al-7) 


N. 


xz 




Is 


Is 


N 


N' 


xz 

5 + 1 
c z 


Is 


1 - 


^'xz^‘ 

I I 
X z 


1 - 


OxP 

I I 
X z 


etc. 


In hovering, the helicopter dynamics are desirable to be decoupled 
into the longitudinal, vertical, lateral and yawing motions, though these 
conditions change in transition and forward flight. The dynamic equations 
to be decoupled are rewritten in state variable form in Eq. Al-8, The broken 
lines make 4x4 submatrices which clarify mam diagonal and off diagonal terms. 
The uncoupled equations with the noise terras added become as given in Eq. 


Al-9. 
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In near hover flight one often assumes u - O^Tg " 0, and the following 
stability derivatives are often neglibly small and often assumed to be 
zero or: 


X = M 
w w 


= Z = 0; 
u 




= Z. 


= 0 


is 


P 






= 0 


With these assumption Equation Al-9 then becomes : 


r n 


— 








1 

- 



» 

u 


X 

-g 

X 

f 

0 

' 0 

0 

0 

1 

0 


u 


u 

q 

i 





1 




« 

0 


0 

0 

1 

1 

1 

1 

0 

‘ 0 

0 

0 

1 

1 

1 

0 


e 

* 


M 

u 

0 

M 

q 

t 

t 

0 

' 0 

0 

0 

1 

1 

0 


q 

— — 


- - - 

- - - - 

- - -* 


« _ - 
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- - - 

- ~ - 

~l ■ 

- - 


" “ 
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= 

0 ' 
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^ ^ . 
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1 Y 

' 'V 
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p 

1 

] 

1 

Y 

r 


V 



0 

0 

0 

1 

1 

0 

' 0 

0 

1 

1 

1 

0 


4) 

• 

p 


' 0 

0 

0 

1 

1 

1 

0 

' L- 
, ^ 

0 

P 

i 

\ 

I 

0 


p 


- 

» _ 

— 



"^1 

- - - 

> 




r 




r 



0 
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-0 

> 0 
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(continued next page) 
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Helicopter Transfer Functions 

For conventional analysis, the equations o£ motions are often written 
in the Laplace domain. The longitudinal and lateral uncoupled equations of 
motion for the control rate derivatives and gust inputs are: 

Longitudinal -Vertical Equations of Motion 
Assumptions 



Lateral-Directional Equations of Motion 
Assumptions 



144 





The basic equations of motion for a helicopter near hover withT =0; 

Uo = 0 

Longitudinal : 
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Lateral : 




The analog computer diagram of the longitudinal and lateral equations of mo- 
tions are shown in Illustration Al-1 and Al-2. 
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Illustration A-la. Lateral Equations 
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Illustration A-lb. Longitudinal Equation 



Input-output Relationships (Hover) 


The uncoupled sets of equations can be used to determine the input- 
output relationships and the effects of gust inputs. 

The longitudinal equations are : 


u 


X 

L+ 

-g 

X ~ 

q 


u 


X 1 
^is 


-X 

u 

0 


0 

0 

1 


9 

+ ! 

0 

Bi - 
Is 

0 

q 


1 

0 

M 

q_ 


q 

1 

M 

B, 

L Is, 


-M 

u 


u_ 


(Ai-13) 


The basic equations are : 


x+Xx-g0 + M6=X„B, -Xu 
u ^ q S^gls u g 

Mx+0 +M9=M„ B, -M u 
u q B^g Is u g 


CM -14) 


or 


s-X 


u 


X^s+g 


~M s -M s 

u q 


u 


X. 


Is 




Is 


Is 


-X 


u 


-M 

u 


u 


g 


(Al-15) 


The characteristic equation is: 


A = s(s-M ) Cs-X ) + M X s + gM 
^ q-^ ^ u’^ u q u 

= s^-OI +X )s^ + CX M +M X )s + gM 
'■qu'^ ^uqug-' 


(Al-16) 


N 


,u 


Is 




M_ s Cs-M J 

Bis q 


CAl-17) 
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n: 


Is 


s-X 

u B, 


-M 


u 


M, 


Is 


(Al-18) 


B 

M„ (s-X + ) 

^is ^B, " 

Is 


The Input-Output Relationships are; 

M X 


u 

Is 


Is 


M, 


Is 


g) 


"B 


Is 


s^-(X +M )s^ + (X M -M X )s + M g 
'• u q*^ ^ u q u q-^ u 


(Al-19) 


Is 


0 

Is 


u 

^is 


M_ (s-X + ^ 
®ls ^ ^B 


Is 


M 3 
u"^ 


Is 


s^-s^(X +M ) + (X M -M X )s + M g 
^ u q^ ^ u q u q u 


(Al-20) 


(Al-21) 


_2_ 

Is 


^B 

Mj, s(s-X + M 3 
B u M u-' 

Is 


s^~s^(X +M 3 + (X N -M X 3s + M g 
^ u q*^ ^ u q u q"^ u® 


(Al-223 



,-I5.Q 



CAl-23) 


Ji, 

U 


g 


U 


g 


„ MX 

-X„(s^ - (M - ^1s -^g) 


u 


Hi. 


M 

i 

Xu 


S^-S^OC +M 3 + rx M -M X 3s + M g 
'■ u q'^ u q u q u* 


. MX M 

J u u 

sfs^-s^CX +M 3 + fX M -M X 3s + M g] 
^ ^uq'^ ^uquq-^ u 


CAl-243 


The lateral equations are; 


r" n 


— 


— 


— - 





• 

V 


Y 

V 

g 

Y 

P 


V 




~Y 

V 

<p 

= 

0 

0 

1 



+ 

0 

A •!• 
Is 

0 


. 

h' 

V 

0 

h' 

P 


P 

_ _ 




-L" 

V 


V 

g 


(Al-253 


The basic equations are; 


V = Yv+g(j) + Yp + Y. A, . 

V A^^ Is 

<^ = Lv+ Lp + L. A^ 

V A^^ Is 


or 


CAl-263 





V 


Is 

A, ■*■ 

-Y 

V 






Is 



S(5-Lj) 


# 


ls_ 


1 

t 

1 


CAl-273 
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The characteristic equation is : 


A = 


s ^ -Cg^YpS^) 

sCs-Lp) 


- (L +Y )s^ + CL Y -L Y )s - L g 

p V P "v V p-^ V® 


N 


Is 


Is 


Is 


■Cg+YpS) 


s Cs-L^) 


2 

Y. [s -CL„. ^YJs + ^ 


Is 


P Y 


A, P 

Is 


g] 


Is 


N 


A, 


Is 


s-Y 

•>r 

Y, 

V 

^Is 

-L 

V 

^A 

^S 

\ 

Is 

[s-Y T 

V 


Is 


L ] 

V-* 


Is 


V 


N 


-Y^ - Cg+YpS]| 


V V 


(AI-28) 


CAI-29) 


CAI-30) 


(Al-31) 


152 



ORIGINAL PAGE IS 
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N.'f 


s-Y -Y 

V V 


-L -L 

V V 


CAl-32) 


= -L (s-Y tY ) = -L s 

V v' V 


The input-output relationships axe: 

Y. - ^YJs + ^g] 


-^“P ■ Y, ‘p- • Y, 


Is 


Is 


Is 


Is 


s[s -CL +YJs^ + CLY-LY)s-LgI 

P V p V V p 


(Al-33) 


V 

Ss 


Y^ [s -(L, 
Is 


Is 


P Y 


Is 


Y„)s + 


'Is 


g] 


Is 


[s^-(L +Y )s^ + (L Y -L Y 3s - L g] 
^ ^pv ^pvvp V ■' 


(Al-34) 


Is 


* r^V 


Is 


Is 


s^ - (L +Y )s^ + (L Y -L Y )s - L g 

^ p v' p V V P^ V* 


(Al-353 


^Is 


Y, 


Is 


L ) 

1 T •• 


L. s(s-Y + ^ 

*ls “ \ “ 

Is 


s^ - (L +YJs^ + (L Y -L Y )s 

^ P V-' ^ P V V P 




(Al-36) 


Y = 


-Y [s‘"-CL - Y )S + ^ 

V ^ p Yv P Y,- ^ 


P Yy p-^ 


L5L. 


s[s -(L +Y )s^ + (L Y -L Y )s - L g] 
p u'^ p V V p*^ ‘V 


(Al-373 
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(Al-38) 


_v 

V 

g 


-v^Cs^-CL fg) 

^ u ^ u 

s^-(L +Y + CL Y -L Y Ts - L g 

^ p V'^ '' p Y V p-^ 


$ 

r 

g 


-L s 
V 


s^-(L +Y,Js^ -h CL Y -L Y )s - L g 

^ P V-^ ^ P V V p-^ V* 


CAI-39) 


V 

g 




s^-CL +Y,Js^ + CLY-LY)s-Lg 
^pv pvvp"^ V* 


CA1-40D 


The vertical equation of motion is : 


w 


= Z w 
w 


+ Z, 


^ - 


z w 

w g 


CAl-41) 


The input-output relationships are: 



CAl-42) 



sCs-Z ] 
'• w' 




CAl-43) 


w _ 

Wg “ s-Z„ 


w 

g 


CAI-44) 



2 


W 

g 


-X 


w 


s fs-X ) 


w 

g 
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CAl-45) 


The yawing equation of motion is : 


r = N r + Np 6^ 
r 


CAl-46) 


The input -output relationships are: 


R 




(Al-47) 




s (s-N^) 


(Al-48) 
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Table A1 . Numerical Data 


In this table the numerical data used in the example is listed. 
The data is comparable to a heavy lift helicopter CH54B gross weight 
47000 lbs. 

The inertia are : 

M = 1460 slugs 

I = 47,020 slug £t^ 

= 192,700 slug ft^ 

I = 164,100 slug £t^ 

I = -19,570 slug £t^ 

xz 

The full control range is: 


6 

c 

(41.0) 

.3 rad 

^s 

(40) 

.28 rad 

^IS 

(28.5 

.455 rad 


(13.6) 

.54 rad w/o.t. with over travel T.R. 


The limitations on the controls are (in degrees) : 
Cyclic 26° 

Lateral 16° 

Collective 6.6° - 23.8° 

Directional 31° C20° without over travel) 


-1-56- 



Table A2. Numerical Values CH54B Hover O.G.E 


The dimensional mmerical values are normalized by the inertia 
G.W. = 47,000 lbs. 



U 

V W 

p4 

q=0 r=iit 

X 

-.0169 

-.0008 .0057 

-2.3082 

.6612 -.2853 

Y 

.0016 

-.0405 -.0160 

-.8945 

-2.1780 1.2315 

Z 

.0044 

-.0150 -.2689 

- . 0244 

.2014 2.6315 

L 

.0012 

-.0153 -.0032 

-1.0937 

-1.2943 .3496 

M 

.0024 

.0003 -.0008 

.3513 

-.3569 -.0072 

N 

-.0002 

.0088 .0028 

-.0147 

.0750 -.5199 


6 

c 

A. 

Is 



X 

6.1041 

-1.2890 

+35.6438 

-.0159 

Y 

-18.0959 

+36.1232 

+1.1884 

23.952 

Z 

-292.260 

-2.2376 

1.0596 

-.0339 

L 

-5.387 

23.4155 

-.1215 

6.5737 

M 

-.1719 

.0648 

-5.6617 

-.1756 

N 

11.627 

.0521 

.2236 

9.5612 
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Table A.3 . Numerical Values CH54B Airspeed 60 Knots 


The dimensional numerical values are normalize.d by the inertia 
G.W. - 47,000 lbs. 



U 

V 

W 

p=(j) 

q=6 

r=ij; 

X 

-.0221 

-.0028 

.0049 

-2.073 

1.515 

-.2431 

Y 

.0027 

-.0717 

-.0286 

-1.7885 

-2.024 

1.990 

Z 

-.100 

-.030 

-.5533 

-1.3915 

1.0951 

2.1848 

L 

.0016 

-.0168 

-.0046 

-1.402 

-1.2729 

.5149 

M 

.0020 

-.0003 

-.0016 

.2814 

-.4254 

.0228 

N 

-.003 

.0144 

-.0002 

.0297 

.2506 

-.8294 


X 

6 

c 

-5.4716 

A, 

Is 

-2.4373 

"is 

33.1214 

-.3398 

Y 

-15.229 

35.353 

4.720 

20.842 

Z 

-304.04 

-1.440 

57.749 

.3203 

L 

-2.2480 

22.82 

.9092 

5.7422 

M 

1.2833 

.2814 

-5.355 

.0710 

N 

7.6112 

-.1179 

.224 

-8.3729 
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APPENDIX B1 


DECOUPLING METHOD (LEAST SQUARES) 

Decoupling Using the Least Squares Method 

\flien the number of outputs of interest exceed that of control in- 
puts, the least squares method is employed by minimizing values of 
coupling terms. Consider the following algebraic matrix equation, 

^ = B CBl-1) 

where 


A : n X ra matrix 
B : n X 1 matrix 
X : m X 1 unknown matrix On<n) 

An approximate solution is obtained by minimizing the following: 

Minimize [ [AX-B [ | CBl -2). 

where | [ • [ | stands for an euclidean norm given by^ 

[[aX-B[[ = {Trace[CAX^BX(AX-Bl'^]>'^ — CBl-3) 


Then, the solution becomes: 

X = ca'^ai'^b'^b 

And 

AX = A(a'^A)~^b'^B = B 
The proof is as follows: 

The solution matrix X is given by: 


CBl-4). 


CBl-S) 


^ [[AX - B[| = 0 


CBl-61 



CBl-7) 


3x 


^ {Trace[CAX-B) (AX-B )^]}2 = 0 


A (AX-B) = 0 


(Bl-8) 


T -IT 

X = (A A) A B, if A has no same column. 


CBL-9) 


This proof uses the following properties of the Trace: 


~ Trace [Ax] = a"^ 


CBl-10) 


3x 


^ Trace [XX^J = 2X 


(Bl-ll) 


A weighting for rows of the matrix equation CB3-1) can be used to in- 
crease the relative accuracy of approximation. The norm becomes: 

1 

1 IqCAX-B) I [ ={Trace[CAX-B)^Q^QCAX-B)J}2 CBl-12) 

where the weighting matrix Q is given by: 


Q = : 


0 0 

22 


0 - 


‘nn 


(Bl-13) 
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The solution is; 


X = (a'^QAX ^ a Q B CBl-14) 

where 

a = Q^Q CBl -15) 

The limitations o£ this method can be illustrated as follmtfs: 

The dynamic equation of the system is: 

X = Ax + Bu C51-16) 

y = X CBl -17) 

(A) If B^Ci-th row) = 0, the corresponding i-th state equation 
should not include coupling terms. 

CB) If y = pXj the output matrix P has only one non-zero element 
in the same column so that desirable state matrix A* and control matrix 


B* can be defined. 



APPENDIX B2 


DECOUPLING OF HELICOPTER DYNAMICS 

Decoupling in the multivariable systems has been studied by many 
authors. Morgan (Ref. Bl) first posed the decoupling problem. Falb 
and Volvoich (B-ef- B2] gave necessary and sufficient condition for de- 
coupling. Gilbert and some authors (Ref. B3, B4, and B53 extended this 
method. However, these methods are still cumbersome in algebraic manipu- 
lation and are often not applicable to control problems where the number 
of controls is less than the number of control system outputs. 

In this work is used a decoupling method of least squares. 

The dynamical equations are given by: 


X = Ax + Bu 

CB2-1] 

y = X 

(B2-2) 

u = Ckx + C6 

(B2-3) 


where 

X : n - state vector 

y : n - output vector 

u : m - control vector 

6 : m - control input vector 

A : n X n system matrix 

B : n X m control matrix 

C : m X m control decoupling matrix 

K : ra X n decoupling state feedback matrix 
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Substituting Eq^ B2-3 into B2-1 and obtain the following equation; 


X = CA + BCKlx + BC6 CB2-4) 

y = X (B2-5) 

Then, the block-diagram for decoupling is shorn in Illustration B2-1. 



State Decoupling 
Matrix 


Illustration B2-1. Block. Diagram for Decoupling 

IVhen the desirable decoupled control matrices are respectively B* 
and A"^, then the following relations are satisfied: 


BC = B* 


(B2- 6) 


A r BCK = A* 


(B2- 75 


From these equations, using a pseudo-inverse matrix which is a solution of 
least squares method, as shora in Appendix Bl, one obtains the decoupling 
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matrices C and K 







T -1 T * 
C = CB B) b' B 


GB2-8i) 


T T -ITT * 

K = CC B BC) e B CA - A ) 


• *T —1 *T 
f CB BD B CA - A 1 


(B2-9) 


Now, tlie decoupling matrices (C,YiX decouple the system completely or 
approximately. The practical example o£ the almost decoupled system is 
shown in Illustration B2-l< This method is used for decoupling of 
the helicopter dynamics at hover and approach phase. 

Decoupling the helicopter Ca-t hover) 

The dynamic equation is reiirritten, using a numerical hover example, 
as follows: (eg. 2.8) with cancelled inertia coupling term) 


u 


-.0169 

-32,2 

.661 1 

.0057 

-.0008 

0 

-2.31 J 

-.285 

0 




1.0 1 


I 


1 


• 

q 


.0024 

0 

-.257 1 

-.0008 

.0003 

0 

.352 * 

1 

-.0072 

w 


.0044 

0 

.201 1 

-.269 

-.015 

• 0 

-.0244 ] 

2.63 

• 

V 


.0016 

0 

-2.18 1 

-.016 

-.0405 

32.2 

-.895 * 
< 

1.231 

i 



0 





1.0 ' 
1 


p 


.0014 

0 

-1.39 ‘ 

- . 0046 

-.0200 

0 

i-H 

• 

»— 1 

.595 

• 

_ r _ 


-.0004 

0 

,241 1 

,0034 

.0112 

0 

.120 ' 
1 

-.591 


__i64 
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35.6 ' 

1 

6.10 

' -*1.29 
1 

' *.016 
1 


1 

0 * 
] 

0 

1 

* 0 
1 

1 

* 0 
1 


1 

-5.66 ' 
1 

-.172 

j 

' .065 

1 

i 

' -.176 

1 

+ 

1.06 ' 

-.292 

i 

' *2.24 

I 

"i ~ ^ “ 

‘ .034 

I 


— - r 

- - - - 

r 



1.19 1 

-18.1 

* 36.1 
1 

* 23.9 
1 


0 * 
V 

0 

* 0 
1 

* 0 
1 


t.223| 

-10.8 

' 24.6 

I 

' 11.1 
1 


_ _ _ 

- - - - 


1- - - - 


.246| 

12.9 

* -2.87 

I 

1 

* -10.9 
1 

1 


And the desirable control matrix B* and the desirable system matrix A* 
are assumed to be given by Eq. and Eq. CB2-12). 


35.6 

0.0 

I 

J 

0.0 , 

0.0 

0 

0 

1 

I 

0 ! 

0 ' 

-5.66 

0.0 

I 

0.0 , 

0.0 

0.0 

-292. 

J 

1 

0.0 , 

0.0 

0.0 

0.0 


36.1 , 

0.0 

0 

0 


0 , 

0 

0.0 

0.0 


24.6 1 

0 

0.0 

0.0 


0.0 , 

-10*9 
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A* 










-.0169 

-32.2 

.661 

0.0 , 

0.0 

0 

0.0 1 

0.0 



o 

o 




1 

1 


.0024 

0 

-.257 

o.a 1 

0,0 

0 

0.0 ' 

0.0 

0.0 

0 

0,0 

-.269 1 

0.0 

0 

0.0 ' 

0.0 

0.0 

0 

0.0 

0.0 1 

-.0405 

32.2 

-.895 1 

0.0 







1.00 1 


0.0 

0 

0.0 

0.0 ' 

-.0200 

0 

-1.14 

0.0 

0.0 

0 

0.0 

0.0 J 

0.0 

0 

0.0 ' 

-.591 
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These matrices B* and A* are chosen in such a way that characteristics 
of the decoupled helicopter are similar to those of the basic helicop- 
ter. After being decoupled, the characteristics would be improved by 
the conventional technique or the optimal control technique. 

After simple calculations following Eq. (B2-S) andCpq. B2-9), 
one obtains the decoupling matrices C and K given by Eq. (j2-l33and 
Eq. CB2-14). 


.998 

-.176 

.042 

-.025 

.004 

1.000 

- .009 

.004 

-.014 

-.159 

1.152 

-.576 

-.005 

1.085 

-.266 

.977 
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-. 168 . 10 “^ 

0.0 

-. 611 - 10 “'* 

1 

. . 155 - 10 “^' . 232 * 10 "“ 

0.0 

-. 648 - 10 “^ 

761 * 10 "^ 

-.ISl'lO"** 

0.0 

. 692 - 10 “^ 

\ 

I-. 237 - 10 “®' .' i ! 4 - 10 “'* 

0.0 

. 834 ‘ 10 "'* 

1 

i -. 900 - 10“2 

. 478 - 10 ’'* 

0.0 

-. 590 - 10 “^ 

.-. 358 - 10 “^ ' . 166 - 10 “^ 

0.0 

. 170 - 10 “® 

r 

i . 308 - 10 “^ 

. 240 - 10 “'* 

0.0 

-. 255 - 10 “^ 

I-. 531 - 10 “^ ‘-. 103 - 10 “^ 
1 

0.0 

. 110 - 10 “^ 

r - • - “ " 
I . 865 - 10 “* 

Lr 
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The member of ths diagonal submatrices of the matrix K are about a factor 
10 ^ less than those of the same column and rows. This means that the 
feedback, gains of the main diagonal sufamatrices, which artificially gov- 
ern characteristics of the helicopter dynamics, have no effect on de- 
coupling with exception of the yawing motion. However, the latter is 
sufficiently stable as is apparent from Eq. (B2-14). This indicates that 
there exist free choice of main diagonal feedback gains without deficiency 
in decoupling. 

Substituting Eq. (B2-14), (B2-13) into Eq. CB2-4), gives the almost 
decoupled helicopter dynamics as follows : 


*u 


-.0169 

-32.2 .661 

1 

I 

.0001 

’ -.0001 
1 

0 

-.0056 

i - 

, -.0151 


u 

8 



1.0 

1 

# 


t 

0 

0 

1 

1 


0 



.0024 

-.266 

1 

.0000 

] .0001 

0 

-.0174 

-.0403 


q 

w 

- 

,0000 

^ 

-.0005 

1 

-.269 

^ .0006 

0 

.0013 

.0030 


w 

V 


-.0001 

-.0203 

j 

j 

-.0021 

'-.0381 

1 

32.2 

-.836 

.137 


V 

4 > 


0 

0 



1 

1 


1.0 



4> 

P 


.0002 

.0279 

J 

.0029 

'-.0233 

JL . _ . _ 

0 

1 

-1.22 

-.189 


p 

r 

• m 


-.0001 

-.0159 

1 

1 

-.0016 

' .0019 
1 

0 

.0466 I 

t 

-.482 


r 
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In order to study the degree o£ decoupling for a control input, an 
analogue simulation was done as shoivn in Figures- C'B2-1] to (B2-8) . The 
decoupling reduces the longitudinal/lateral coupling terms by a factor of 
10 or more. The coupling between the vertical mode and other modes also 
decreases, showing significant reduction between the vertical and yawing 
mode. The yawing response is improved for the control input, while its 
coupling with the other axises decreases significantly. A change in some 
control derivatives signs accompany turning direction of initial motion 
opposite to that of basic derivatives . An evaluation of this sign change 
may require a simulator test by a pilot in the case of an unstable helicopter. 


Decoupling the Helicopter Cat approach phase) 

The dynamic equation, using the numerical example for the approach plane, 
becomes as follows: (with cancelled inertia coupling term) 


• 

u 


-.0221 

-32,2 

1.S15 ' 

t 

.0049 

-.002 

0 

-2.73 

,-.2481 

• 

9 


0 

0 

1.0 ‘ 

0 

0 

0 

0 

1 

( 0 

• 

_q_ 


.0020 

0 

-.4251“ 

.0016 

.0003 

0 

.281 

I-.228 

• 

w 


-.100 

_ 0 _ 

7.845 ' 

-.5333 

-.030 _ 

0 

-1.3915 

,2.1898 

• 

V 


.0027 

0 

-2.024 ' 

f 

-.0286 

-.0711 

32 

,2 -1..788 

,1.990 



0 

0 

0 ' 
1 

0 

0 

0 

1.0 

1 

1 0 

• 

p 


.0629 

0 

-1.449' 

-.0047 

-.0113 

0 

-1.488 

',2.169 

• 

r 


-.0003 

0 

1 

.423' 

.0057 '' 

. 0172 

0 

.207 

I-.937 _ 
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33.12 

-5.470 

-2.44 

-.361 

0 

0 

0 

0 

-5.35 

1.28 

.28 

.076 

57.75 

-304.04 

-1.44 

.32 

4.72 

-15.221 

35.35 

20.84 

0 

0 

0 

0 

.86 

-5.69 

24.06 

9.70 

,122 

8.29 

00 

-9.53 
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and the desirable control matrix B* and the desirable system matrix 


B* = 


1 given by 

Eq. 

(B2-18 ) 

and Eq. 

33.12 

0.0 

0.0 

0.0 

0 

0 

0 

0 

-5.35 

0.0 

0.0 

0.0 

0.0 

•304. 

4 0.0 

0.0 

0.0 

0.0 

35.35 

0.0 

0 

0 

0 

0 

0.0 

0.0 

24.06 

0 

0.0 

0.0 

0.0 

-9.53 
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.0221 

-32.2 

1,515 

1 1 
1 0-0 1 

0.0 

0 

0.0 

1 

1 0.0 



1.0 

1 1 

I ! 




1 

1 

.0020 

0 
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i ] 

1 0-0 1 
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0 

0.0 

1 

1 0.0 
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0 
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p 

1 -.533^ 
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0 
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1 0.0 
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— - 

- r 1 

- - - - 

- - - 
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r ~ - 

0 
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1 0-0 1 
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32.2 

-1.7SS 

1 0.0 




1 1 
1 1 



1.0 

1 

1 

0.0 

0 

_ 0.0 _ 

1 1 

_l L 

-.0113 _ 

_0_ _ 

-1.488 

1 

.i 

0.0 

0 

0.0 

i ! 

1 0-0 ! 

0.0 

0 

0.0 

1 

1 -.937 


(B2-18) 

These matrices and A* are chosen in such a way that characteri'stics 
of the decoupled helicopter are similar to those of the basic helicopter. 
After being decoupled, the characteristics would be improved by the con- 
ventional technique or the optimal control technique. 

After simple calculations following Eq. (B2-16) and Eq. (B2-17), one 
obtains the decoupling matrices C and K given by Eq. CB2-21) and-Eq. (62-22). 


C = 


1.02 

.17 

.08 

-.034 ! 


,19 

1.03 

.010 

.003 


-.1 

-.14 

1.15 

-.5 

(B2-19) 

.13 

.95 

-.30 

.986 
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0.0 

0.0 

0.0 I . 151 ^ 10 "^' 

- r - 

. 102 - 10 "** 

0.0 

-. 816 - 10 “^ 

.-. 621 * 10 "^ 

_. 328 . 10 ~® 

0.0 

i 

-. 360 . 10 "^ . 0 ' 

.986 ‘ lO "** 

0.0 

. 457 - lO " 

.-. 717 - 10 “^ 

. 903 - lO "** 

0.0 

. ^ ^ ^ r 

-. 581 - 10 "^ i ~. 603 - 10 "^ ] 

0.0 

0.0 

0.0 

' . 670'10 ^ 

. 314 - 10 """ 

0.0 

-. 443 - lO "^ t -. 388 ' 10 “® 

. 181 - 10 '^ 

0.0 

. 220*10 ^ 

> 0.0 
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Substituting Equation B2-9 and B2-20 into Equation B2-5, gives the almost 
decoupled helicopter dynamics as follows. 

u 

e 
q 

w; 

V 

4 > 
P 
r 

CB2-21 ) 
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0 
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0 
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-.0109 

0 
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0 

0 
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0 

0 

0 
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33.02 

.077 

.041 

-.210 

0 

0 

0 

0 

-5.24 

.43 

-.137 

.106 

■>0008 

-304.04 

.0003 

-.0018 

.01 

.077 

34.76 

2.10 

0 

0 

0 

0 

1,39 

.11 

24.8 

-.029 

.79 

-.057 

-.459 

-7.88 




6 
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Note that in the "decoupled" equations, the coupling terms are an order of 
ten smaller, compared to those of basic equations. 

In the numerical values of the above equations is used with 

U = 0. If Z -U is used with U = 60 knots/hr, the feedback gains k. 

O R. O 0 IX 

remain the same with the exception of which becomes equal to -100 i.e., 
directly proportional to speed. 
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Evaluation of Decoupling 


The root configuration for basic, decoupled and uncoupled simpli- 
fied dynamics are shown in Figure 2. The roots of the decoupled longi- 
tudinal equations are very close to those of the uncoupled equations. 

However, decoupling of the lateral equations shifts one of the modes to 
the right hand side in the s plane, resulting in an unstable mode. 

In order to study the degree of decoupling for a control input, an 
analog simulation was done in Reference (32) . The decoupling reduces con- 
siderably the longitudinal/ lateral coupling terms. The coupling between the 
vertical mode and other modes also decreases, showing significant reduction 
between the vertical and yawing mode . The yawing response is iraporved for 
the control input, while its coupling with the other axes decreases signifi- 
cantly. A change in some control derivatives signs accompany turning direction 
of initial motion opposite to that of basic derivatives . An evaluation of 
this sign change may require a simulator test by a pilot in the case of an 
unstable helicopter. 

Consideration for Practical Design 

In general, the stability derivatives and control derivatives are a 
function of air velocity, altitude and time. Therefore, each element of the 
decoupling matrices [C,K] must be given by: 



= C (V,h,t) 
-*• J 



K.^(V,h,t3 


(B2-23) 
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However j the derivatives are determined at a given air velocity and 
altitude. In practice, the altitude and time dependence can often 
be neglected and the derivatives are only a function of air velocity 
or: 



= C. . 
^3 


CV) 


K.. 

13 


= K. 


13 


Cv^ 
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IVhen the derivatives between different velocities can be approxi- 
mated by straight lines or simple curves variations in C. . and K. . 

13 13 

would be easily programmed. 

The decoupling matrices C and K are constructed as shoAvn in Jl lustration B2- 
3. The assumptions for the design of the helicopter controller are 
as follows: 

Cl) The state variables u, v, w, 9, q, 4'^ P aud r can be mea- 
sured by sensors or approximated. 

Cni The fly-by wire system would be desirable. 

CIII) The measurement of air velocity V is possible. 
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APPENDIX Cl 


SENSITIVITY ANALYSIS 

Sensitivity of CRaracteristic Roots to Variation. 'in the System 
Coefficients 

THe system equation vd.th, control free is given by (pl-1)- 
x(t) = Ax(t) 

where 

; n-state vector 
A ; n X n system matrix 

the n-th order characteristic equation is given by 

isi - A| =0 Cci-2] 

The sensitivity is defined as 


where. 


k H 

S.. = a. . 

1 ] ij ia. . 


CCl-3) 


s^. 


a. . 
Aa. . 

Ar, 


sensitivity of k-th characteristic root 
to a. . variation 
non-zero element of A 

variation in a. 

k-th root of the characterxstic equation 


And, S. , is a vector since Ar, is represented by: 

X 


, Im[Ar, ] 

Ar^ = llAxJ I tan- 


CCl-4) 


Also, the following equations are easily proved with respect to , 
since Trace [A] is the sum of all roots. 


-1-76 



i = j, i,j= 1,2 


(C 1-5) 


n 

T. 

k=l 



a. . 

ij 


n 


S Sj = 0 i j, i,j = 1, 2, n CCl-6) 

k=l ^ 


Furtfiermore, the poles and zeros of the root locus of a. - varia- 

X j 

tion are given by~ roots of (^1-7) and '(Cl-S) respectxvel/. 


si - a| = 0 


cofactor 


C-a. , ) = 

ir 


CCl-7) 

CCl-8) 


These are derived as follows: 

Letting the variation of the system matrix, AA, the characteristic 
equation is rewritten by: 

jsl - A - AA] =5 0 ©1-9) 


This equation can be equivalent to (Cl-10) 
|l - (si - A)"^AA| = 0 


(Cl-10) 


adj CsI-A)AA j _ » 
detCsI-A) ' 


(Cl-11) 
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co£(s-a^^), cofC-a^^), cofC-a^^^), 

cof C-a, , co£ Cs-a„«) , , co£ (-a. ) , 


, co£C-a^^) 

, co£C-a^2) 


co£C-a^^), cofC-a^^) 


co£C~a 0 
^ nj 


cof(-aj^), cofC-a^jj), 




l£ one assumes that 


then: 


adj CsI-A)AA = 


0 0. 


0 


. . .Aa^j.co£C-a^j^) 
Aaj..cofC-aj2) 


,4ay.cofC-ay) 


0 Aa^j, co£C-a^^). 


• adjCsI-A)AA 

’ - det[sI-Aj 


= 1 - 


Aa^^co£C-a^^) 

det[sI-A] 


That is, the characteristic equation is: 

Aa. . co£ C-a. . ) 

1 13 13 _ n 

“ det[sI-AJ 


The poles and zeros are given from p 1-16) 
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Sensitivltes of the Characteristi^c Roots to the Stability Derivatives 


TPie contributions o£ stability derivatives to control and dynamic 
stability are shown in Figures P-l}_ to (t-361j with, mostly expanded 
sensitivities. By these figures w;e can appreciate the following: 

Cll IVhich roots Cor modest axe most sensitive to the specific 
stability derivative? 

C21 IVhich roots Cor modes) are more sensitive to which stability 
derivatives? 

C3) Existence of zeroes on real axis can be knoivn. 

C4-1 From Cl) C2) we can infer degree of coupling. 

With respect to L and N .root locus are shown with reducing L and N 
^ r u ° r u 

in Figure (C-37) and Figure (C-oS), in order to study reduction in 

stability due to decoupling. 
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APPENDIX D1 

review of linear optimal regulator 


Optimal Time -Invariant Linear Regulator (Deterministic) 


The general results of optimal time-invariant linear regulator 
problems are summarized in this section. 

System Dynamics 

Given is a deterministic, time-invariant, linear dynamical system 
with n-dimensional state vector (or variable) x(t) and m-dimensional 
control vector (ot variable) u(t) . 

xCt) = AxCt) + Bu(t) (Dl-1) 

xCt^) = 0)1-2) 

where A and B are n x n and n x m constant coefficient matrices, and x 

o 

is the initial state vector at time t^. Usually, the dimensional of state 
vector is larger than that of control vector (n>m) , 

The p-dimensional output vector yCt) consists of a linear combina- 
tion of the components of state vector xCt) : 

y(t) = Hx(t) (Dl-3) 

Without loss of generality, the p x n' output matrix H has its maximum 
rank p Cp^^) . The output vector yCt) can be regarded as p independent 
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measurements o£ the system state x(t). Note that, in many or most o£ 
the practical cases, the output matrix H can be an n x n identity matrix. 
In Illustration (Dl-1) is shown the block diagram o£ the system. 

Throughout this work, the lower case character usually stands for 
a column vector and the upper case for a matrix. The argument ”t” usu- 
ally means time and C‘) represents the time derivative of the argument. 

Cost Function 

By properly chosen output weighting matrix Q and control weighting 
matrix R, the quadratic cost function without cross coupled term of the 
output vector and the control vector is expressed as follows 

J(x^, u(t)) = y* [y'^(t)OyCt) + u^Ct)RuCt)]dt 

^0 CD 1-4;) 

where the matrices Q and R are pxp and mxra constant symmetric matrices, 

T 

and throughout this report, the superscript T CC') ) stands for the 
transpose of vector or matrix. 

The pxp symmetric output weighting matrix Q and the mxm symmetric 
control weighting matrix R are both assumed to be positive definite, or 

Q > 0 

} CDl -5) 

R > 0 

Hence, the integrand of quadratic cost function J(x^,uCt3) in Eq, (Dl-4) 
is non-negative for any output y(X) and any control uCt) during the pro- 
cess interval [t ,«). 

0 
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C ontro liability 


An optimal control law u'*' (t) has to be determined, that minimizes 
the quadratic cost function u(tX) sub.j.ect to the system dynamrcs 

represented by Eqs. (pl-1), (Pl-2) and (pl-3). 

The assumption of the positive definiteness for Q and R weighting 
matrices guarantees the convexity of quadratic cost function ti{t)), 

but it does not imply that the minimum of cost function J fx ] is finite. 
The cost function is finite by the "complete controllability condition". 

The condition of "Complete Controllability" for stationary system 
(Eqs. (Dl-1) and (Pl-2) considered here is that the following n rows by 
nxnm composite matrix has the maximum rank (=n) : 


Rank [B : AB 


A^B 


a""^B] = n (Dl-6) 


A matrix pair of A and B satisfying Eq. (Dl-6) is usually said as "Com- 
pletely Controllable pair [A, Bj". 

Complete controllability guarantees the finiteness of the minimum 
cost function uCt)). Also, since the quadratic cost function Eq, 

Cd 1-4). is a convex function, the uCtDD exists and is unique. 

Optimal Control Law 

The optimal control problem can be solved by the aid of calculus 
of variations, or Pontryagin's maximum principle and/or Bellman's dynamic 
programming. Since there is no active constraint involved except the dy- 
namical equation, the simplest straightforward way may be to use the cal- 
culus of variations. The result of this calculation is that matrix K must 
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satisfy the "Algebraic Matrix Riccati Equation". 


KA+A^K - KDK + H^QH = 0 (Pl-7) 

(note H^QH >03 

where D is a n x n constant, synmetric, positive semidefinite matrix 
defined by 


D S C^O) (Dl-8) 

Let the matrix K be a proper solution of Eq, (Dl-7 Then the optimal 
control law is expressed in the linear feedback form. Using Eq. (Dl-15) 
and (Dl-nh 


u-^Ct) = -R"^B'^Kx*Ct) 

= -Gx-^Ct) 

where G is the feedback gain matrix (m x n) defined by 


CD 1-9) 


G = R’^b'^K (Dl-10) 

By using the optimal feedback control law to the plant dynamics Eq. (Dl-1), 
the optimal closed loop system dynamics become 


xCt) = (A-BR'^b'^K) xCt) = FxCt) 


CDl-11) 
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where the asterisks on x(t) have been dropped off, and the n x n matrix F 
is the closed loop system matrix defined by 

F = A-BR'^b'^K = A-DK = A-BG Cd1-12) 

The eigenvalues of matrix F indicate the poles of the optimal closed 
loop system. The stability of the system depends on the signs of the real 
parts of the eigenvalues. There exist many real symmetric solution ma- 
trices CK's) of the algebraic matrix Riccati equation C^q. d 1-7) for given 
positive definite matrices Q and R. It has been shown by R.E. Kalman in 
Ref. (aJ) that if [A, B] is a completely controllable pair, there exists 
a unique solution matrix K of Eq. (Dl-19), which is a n x n positive defin- 
ite symmetric matrix and makes the real parts of all the eigenvalues of F 
be negative (it means a stable system); namely, 

K > 0, unique and symmetric 
for Q>0 and R>0 

and 

R^{X^[*]} means the real part of the i- th eigenvalue of a matrix. 
Therefore, the closed loop system expressed by Eq. (pl-ll) is asymptotically 
stable. Thus, the optimal control law u*Ct) in Eq. (pl-9) exists and is 
unique for given weighting matrices Q and R. 

It is important to note that the matrix Riccati equation (Eq. (P 1-7) 
is independent of the system state and time, and so is the optimal feedback 
gain matrix G in Eq. C^l-10). 
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observed vector) yCt)^ despite the attempt to use the output y(t) as the 
feedback signal. Hence, the optimal controller requires all the infor- 
mation about the system state x(t), which has to be computable out of the 
observed vector yCt). If the output vector y(t) could not generate the 
(unique) state vector x(t) or if the state vector x(t) could not be di- 
rectly observed (qt measured) , then the obtained optimal system would not 
be realizable. 

A "state x^ is observable at t^" if, given any control u(t), there 
is a time such that knowledge of u(t) and the output y(t) over 

is sufficient to determine x^. If the argument is true for every state 
-X at time t , then the system is said to be "observable at t ". If every 
x^ is observable at every time in the interval of definition of the sys- 
tem, then the system is said to be "completely observable". Since the sys- 
tem under consideration here is stationary (constant coefficient), the sys- 
tem (Eqs. (Dl-1), (Dl-2), (Dl-3)) is assumed to be completely observable. 

The condition of complete observability for time-invariant linear sys- 
tem is that the following n-rows by np-columns composite matrix has the 
maximum rank (>n) 


Rank [H 




a‘h":(a'^)^h'^ 


,.T.n-l„TT 

CA ) H ] = n 


(Dl-12) 


A matrix pair of A and H satisfying Eq. C^l-30) is said "Completely Ob- 
servable pair [A,HJ", 
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Minimum of Quadratic Cost Function 


By using the unique optimal feedback control law CEq* the 

minimum of quadratic cost function unfquely evaluated in 

terras of K matrix and the initial state x^; and equals: 


J . (x ; Q, R) = Kx (D l-B ) 

min^ o ’ o 0 

Block Diagram Optimal Regulator 

The block diagram of the optimal system described above is shown as 
Illustration Dl-1. 


Controller 


System 




Illustration Dl-1. Block Diagram of Linear Optimal Control System 
Observability 

Eq. C)l-21) and the block diagram shows that the optimal control 
u*(t) is a function of the state vector xCt]) rather than a function of 
the state vector xCt) rather than a function of the output vector Cor the 
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Note that, combining Eq. (Dl-6) v;ith Eq. (Dl-12), the system is 
assumed to be both completely controllable and completely observable. 

The observability condition Eq, (Dl-12) guarantees that one can 
compute the unique state x(t) from knowledge of the output y(t) and the 
control u(t) over [t^,t]. When the output matrix H is the n x n iden- 
tity matrix (=1^^) > the problem is called "State-Regulator Problem", where- 
as the problem with a nonidentity H matrix is called "Output -Regulator 
Problem" . 

Summary of the Deterministic Optimal Regulator 
The deterministic optimal regulator for a completely controllable 
and observable system can be summarized as follows 

System dynamics : xCt) = AxCt) + BuCt); xCt^] = x^ 

yCt] = HxCt) 

where x(t^ : n dimensional state vector 

u(t) : m dimensional control vector (m^) 

yCt) : p dimensional output vector 

A,B and H: nxn, nxm and pxn constant matrices 

Cost function; = Min f ^ly"^ (t)QyCt) + u^Ct)RuCt)]dt (3) 

min 0 

o 

= x' K X C4 and 11) 

o 0 

where Q and R: positive definite symmetric constant 


matrices 



Variables and equations; 


Number of Variables 

Number of Equations 

x„Ct) : 

n 

Cl] : n 

u(t) : 

m 

CD : P 

yCt^ : 

P 

(3) : 1 

JC-) : 

1 

Cll] or C4) : m 

Q : 


CD : 

R : 

m^ 


K : 



Total: Cn+ra+p+l+p^+m^+n^) unknowns. 

Total: Cn-i-p+l+m+n^) eqs. 

Conclusion: 

1. If Q, R are given. 

then number of variables = ntimber 


of equations. 

2. If Q and R are not 

given, then the number of equations 


is less than the number of unknowns. 
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APPENDIX D2 


OPTIMAL, LINEAR CONTROLLER DESIGN APPLIED TO A SECOND ORDER SYSTEM 
The system equation is represented by the following differential 
equation 


= axCt) + bu(t) CD2-1) 

where 


x(t) = position 
x(t) = rate of position 
uft) = control input 
a = rate damping 
b = control derivatives 
In the state variable form 


xft] = Ax(t) + Bu(t) 

where 


xCt] = 


A = 



B 



uCt) = 5 


CD2-23 


(D2-3) 


(D2-4) 

(D2-5) 


CD2-6) 
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The Optimal Linear Regulator Approach 


The quadratic performance index be minimized is: 


J = 



Cx(t)'^QxCt) + ru(t)^3dt 


CD2-73 


where 

Q : 2 by 2 non-negative, symmetrix weighting matrix for 

qil^ ^12^ ^21 ^ ° states 

r : a scalar weighting factor for control 


The optimal control which minimize (D2-7) is 

C 

u = -Gx(t3 

= - I B'^Kx(t) CD2-S) 

where G is a feedback gain matrix and K is a solution matrix of the 
following Riccati matrix equation. 

KA + A^k - KBB^ + Q = 0 (D2-9) 


Expanding the latter equation gives: 


■•12 


kai - qii 


0 
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Hlcj* 


0 


CD2-10) 


k- k 
12 22 


Cak + k D - q 
12 11 12 


— k k -Cak +k)-q =0 

r 12 22 21 11 21 


^ k ^ 
-r- ^22 


2ak22 - (k-12 + kai ) - q22 = 0 


From the symmetric property o£ weighting matrix Q follows qia = qai. The 

« 

solution of Eq. D2-10 becomes 


Cix - 1 1 r ( C^) ^ F ^ ^ 


qiir 


CD2-11) 


+ ^ 2b/5H + ^ "^223 


Substitution into (D2-S) gives the optimal control 


^ [S X, + G.x,] 
r '■ X 1 X 2-' 


(D2-12) 


where is a position feedback gain and G^ is its rate feedback gain 
given as follows: 


== _ /iH 


(D2-13) 





b r 


22 


CD2-14) 


Note that the weighting factors qi 2 or qzi are irrelevant to the op- 
timal feedback gains. 

From (D2-23 , lt)Z-12), CD2 -13) and (1)2-143, the equation for the 
closed loop system becomes; 


/ 


xCt) + /a^+2b/^ + X + b/Iil 


X(t3 = 0 (D2-15) 


(Vhen this closed loop system is described by 


x(t) + 2?(0^x(t) + 0)^ x(t) = 0 


then 


(D2-16) 


0) 

n 


= U 


j_i 

r 


— '/l 
2 


C = T 2 + 


/ 


qii ^ rqii 


Cqii f 0) CD2-17) 


Closed Loop System Properties as a Function of Control Weighting 
When the control weighting R approaches infinity, CR-x») 


G =0 

X ^ 


= 0 


CD2-18) 
(;D2-19) 

The feedback gains do not all go to zero; therefore the closed loop system 


2a 

b 


(a<o) 

Ca>o) 
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approaches the stable original system or the adjoint system in the case 
of an unstable original system because: 

xCt}' + |a| x(t) = 0 (D2-20) 


always guarantees a stable system. 

Itfhen the weighting for control approaches zero, = 0) , which means 
there is no limitation on the control, the closed loop system approaches 
1132-21) which is called model control equation in (D2-l) . 

xCt) + x(t) = 0 (D2-21) 

V qii 

when the weighting factor q 22 = 0, the undamped natural frequency and 
damping become, from (D2-7) : 


n 


(D2-22) 


C = .707 


which are the second order Butterworth poles , 


The Minimum Value of the Performance Index 


The minimum value of the performance index, in general, is obtained 
by substitution of the optimal control (D2 -8) into CK2-7), This yields: 


1 T 

J . = X Kx 

min 2 0 o 


21L 


CD2-23) 



This is obtained by substituting D2-11 into the latter equation: 


J . = 

ram 


1 

2 



+ 





CD2-24) 


£ 

r 


<122] C 


X20 

xio 


Cxio r 0) 


The formula, in general, is given as a function of initial states and 

solution of the Riccati equation. Minimizing thrs value is done by de- 

/ 

creasing r and also by decreasing qn and qa 2 . Therefore, without any 
constraints, minimi zing J itself is meaningless. Hence, constraints are 
necessary, which leads to the discussion of the OWEM method. CAppendix E) . 
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APPENDIX El 


REVIEW OF OWEM DESIGN i^ETHOD 


The Optimal Weighting Matrices Controller Design Method (OWEM) 


This appendix is a summary of the OWEM design method as presented 
in Refs. (£1 ) and C ^2 ). The results of the method are given without 
proofs; for the details of mathematical derivations and proofs, see Ref. 
C El). 

The goal of the linear optimal control theory is to find a linear 
control law such that the following quadratic cost function is minimized 


J(uCt); Q,R) = Lim 

T->co 


1 

2T 



{y'^Ct)QyCt) + u'^Ct)RuCt)Jdt 

CEl-1) 


subject to :kCt') = AxCt) + BuCt) + CnCt) 

yCt) = HxCt) 

£{nCt)} = 0, £{nCt), n^Cr)} = NdCt-x) 


where xCt) 
uCt) 
yCt) 
nCt) 


n dimensional state vector 

m dimensional control vector 

p dimensional output vector 

£ dimensional stationary, Gaussian, white 
noise vector 


Q : pxp output weighting matrix, weightings for the 
output errors, (positive definite, symmetric, 
constant matrix) 

N : nxn noise intensity matrix 
<S(*) : Dirac's delta function 
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R : mxm control weighting matrix (positive definite, 
symmetric, constant matrix) 

T 

(•) : transpose of (•) 

0{*} : expectation or average of {•} 

The optimal control law u*(t) is obtained in the state feedback form 
(a linear function of xCt)) as follows 


u*(t) = -R"^B^Kx(t) = -Gx(t) (EI-2) 

where (•) ^ represents the inverse of (•), and nxn matrix K is the unique 
positive definite solution of the following so-called "Algebraic Matrix 
Riccati (AMR)" equation: 


KA + a"^K - KBR"^b'^K + h’^QH = 0 (El-3) 

Thus, if the matrices A, B, H, 0, and R are given, the controller can be 
designed by solving Eq. (^1-3 ) ^ a.nd using it in Eq.'CE'i^)- 

Note that the optimal controller can be designed independently of the 
noise intensity (covariance) and the matrix C as long as the noise n(t) 
is stationary, gaussian, white and of zero-mean. 

Denote the corresponding state and output vectors to the optimal con- 
trol u*(t) as x*(t) and y*Ct) respectively. Namely, 


x*Ct) = Ax*(t) + Bu*(t) + Cn(t) (El-4) 

y*Ct) = Hx*(t) (EI-5) 

The equation should minimize J(u(t);Q,R). 
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Assuming the stationarity of the system and the ergodic property, 
the minimum cost JCu*Ct),Q,R) rs obtained as follows (pef. El ) 

J^inCQ,R) = JCu*Ct);Q,R) = £{y*Ct)Qy*Ct) + u*‘^Ct3Runt)3 

= Tr[QY"J + Tr[RU*] CEl-6) 

= Tr[CNc'^K] (El-7 3 

where Tr[*] = Trace of [•] 

= the sum of the diagonals of [•] 

x 

Y^Ct) =£'Cy*C't)y* (t)} : output covariance matrix 
U*(t) = £{u* C_t)u*'^Ct) ^ • control covariance matrix 

The matrices A, B, C and H are inherent to the given controlled 

system and normally available a priori. In order to design the controller, 

the immediate question to be asked is how to choose Q and R matrices. They 

affect the solution matrix K in AMR equation and thereby the feedback gain 
-1 T 

matrix GC=R B K). Hence, the closed loop characteristics and the minimum 
quadratic cost have considerable effects of Q and R. 

First will be investigated how the minimum quadratic cost 
behaves with respect to Q and R: 

Suppose there is an ordered sequence of trial weighting matrices O's. 


[Q^J = 




■ Q., Q. 
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^i +1 “ ^ 0 or positive semidefinite, then [Q^] is called a monotoni- 

cally increasing definite sequence (MIDS) . For an example, consider one 
Q^>o and a scalar sequence [k^J = 1, 2, Zj i, i+1, . Then 

Q. = k.Q 
1 1^0 

[Q^] = Qi, Q2, Q3, •" Q^, •••• 

= Qq, 2Q^, 3Q^, ••• iQ^,CUlDQ^ .... 
is MIDS. Because 

^1+1 ■ ^i " = for all i = 1, 2, 3, ••••. 

Likewise, [R^J and [K^] are defined. 

Using these special sequences, the parametric behavior of the minimum 
quadratic cost summarized below CRef. [El]3. 

A solution sequence of AllR equation corresponds to the above 

defined sequences [Q^J and [R^] . 

Corresponding to the above defined sequence [Q^j , -and [R^] respec- 
tively, one obtains solution sequences of AMR equation. The solu- 

tion sequences {K^}'s can be shown to be also monotonically increasing 
definite sequences [K^J's. This is shown in Illustration E-1. 
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Illustration E-1. Behavior of [K^J with respect to [0^] and [R^] 


The with respect to a monotonically increasing definite 

sequence [K^] also increases monotonically. 



Illustration E-2. Behavior of [K^] (N;Q,R) with respect to [K^] 

The arrow in- the illustration indicates the direction 
that is desirable for a design. 


Define J . (K. ) = Tr [CNc’^K. ] 
min 1 


(J . CN;Q. ,R) 
* mm ^1 ^ 

or 

J CN;Q,R.) 

min^ i'^ 


CEI-8) 


Then if K. - K. > 0, then J . (K. J > J . (K.) 

1+1 1 — ' min 1+1 — min 


(El-9) 
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The behavior of with respect to [Q^] and [R^^J is shorn 
in Illustration E-3. The sequences [Q^] and [R^^J are monotonically in- 
creasing definite sequences. 




Illustration E-3. J . with Respect to [Q.l and [R.] 

rain ^ '■ i-* 

The arrows in the illustration indicate the direction that is desirable in 
a design. 

The raonotonic property of JCN;Q,R) indicates that will not remain 

finite unless the weighting matrices Q and R are limited within some finite 
allowable sets during their selection. V/hen the elemental matrices of the 
allowable sets for choosing the Q and R matrices can be ordered and a 
definite differential matrix among any two of them is obtained; then the 
smallest weighting matrix yields the smallest J . (N;Q,R). Hence, the 
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alloi^able sets for the weighting matrices Q and R should not contain any 
matrices which can be ordered in the sense of definite differential matrix. 
To bound J (N:Q.R) from both below and above, consider Eq. ) 

J^inCN;Q,R) = Tr[QY*] ^ Tr[RU*]. 

The general form to be bounded is 

Tr[A B] where A corresponds to Q or R 

B corresponds to or U* . 

This trace indicates a sum of the products: 

TiCa'^B] = J f (El-iO) 

where A = {a. B = {b. .}. 

ij ij 

For* a sum of productes of finite values, the lower bound is given by the 
arithmetic geometric average, the upperbound by the Cauchy-Schwarz theorem. 
For an example, 

2/aia2bib2 ^ aibi + a 2 b 25 . + a| + b| (El-U) 

where a,b,>0 and aab 2 > 0 for the left hand side inequality. 

In matrix form this becomes [Ref. El -3]: 

IVhen A and B are non singular Cnxn) matrices then 

n |A|'/^-iB|^/^ < Tr[AB]< liAll • l|B|i (El-12) 
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where 


I • I : Determinant 

I j • I I : Matrix norm, 1 |Aj | = Jlr [A^AJ 


The equalities occur at 

~T ~-l 

the left hand side when A = 

the right hand side when A = ^2^^* 

Hence J . (N;Q,R) is bounded by 


1 llQli -MY-jl - ||R|1- ||U*i| 


Eq. E-12 indicates that: 

1. IVhen iQj and |r| are prescribed, the minimum of J . CN;Q,R) with 

min 

respect to Q and R can be obtained with Q = k Y*”^ and R = k.U*~^. 

O ^ 

2. When ||Q|j and |[R(| are prescribed, the maximum of 

with respect to Q and R can be obtained with 0 = and R = ^5^^* 

where k^, k^, k^ and k^ are appropriate positive scalars for Q and 
r to satisfy the prescribed determinants and norms. 


The minimum and maximum values of J . CNrO,R) occur when: 

min^ > .> f 

Min J . CN;Q,R1 when Q = k_Y*"^ and R - k.U*“^ 
min ^3 4 

Max J . CN;Q,R) when Q = k_Y* and R = k,U* 
nixii o o 
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The design method with Optimal Weighting Matrices COWEM) is basi- 
cally a linear optimal quadratic control system design method with calcu- 
lated weighting factors such that the effect of disturbances on the per- 
formance index is minimized. 

The OWEM design method can be summarized as follows: 

The system is assumed to be completely controllable and observable and is 
shown in Illustration (El-l). 



Illustration El-i . Block Diagram of an Optimal Linear Regulator 


System dynamics: 


5c (,t^ = Ax(t) 
yCt) = HxCt) 
where x(t) : 
uCt] : 
y(t) : 

nCt) : 


+ Bu(t) + CnCt) CE1-I4) 

CEI-15] 

n dimensional state vector 
m dimensional control vector Cra<n) 
p dimensional output vector (pfn) 
a dimensional noise vector (ii<n) 
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£{nCt)} = 0, Cov{iiCt)} = N6(t-T) 

A,B,C and H; nxn, nxm^ nx5. and pxn constant 


given matrices. 


CEl-16) 


Cost function : J . (N;Q,RD = Min {Lim ~ [y^Ct)Qy(t)+u"^Ct)RuCt) ]dt} 

uCt) T- -T 

= Tr[CNC^K] 

where Q and R: positive definite symmetric matrices. 

-1 T 

Optimal control law : u’lt) = -R B Kx(t) (El- 19) 

where K: nxn unique positive definite solution of (El- 20) 


KA + a"^K - KBR"^b'^K + h'^QH = 0 
T T 

Covariance, matrices: FX + XF + CNC = 0 (Ref. E2) 


where X = 6{x(t)x'^(t)}, F = A - BR'^b'^K 

Y = £{y(t)y'^Ct)} = HXh"^ 


u = £,{u(t)u”^(t)} = r"^b'^kxkbr"^ 


(El- 20) 
(El-21) 

(El-22) 

(El-23) 
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In the above equations, the number of equations is p +m less than 
the number of unkno^ms if Q and R are not given. To augment the lacking 
number of equations, an atnciliary performance index is minimized with re- 
spect to Q and R. 

Auxiliary performance index; 


Min J . (N;Q,R) = Min Tr[CNC K] 

Q,R Q,R 

subject to KA + a’^'K - KBR'^b'^K + h'^QH = 0 


(El- 2 4) 
(El- 2 5) 
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(£1-26:) 


|r| 

Constraints ou Q and R ; jQj = 1, |q|= [j^j = p 2 

Optimal weighting matrices Q and R : 

Q = ^IhPh'^ j • (El-27) 



b'^kpkb 


I b'^kpkb I 



(El- 28) 


where FP + PF^ + CNC^ = 0 P E X in this case. (£1-^) 

F = A - BR"^ b'^K 

Conclusion ; 

If is given, then Q and R can be determined provided the noise 
intensity matrix N is known. If the noise intensity matrix (or covariance 
matrix) N cannot be estimated, then the approach in this summary cannot be 
used. 

It has been sho™ in Ref. (El ) that the introduction of a new per- 
formance index independent of N can remove this deficiency. This performance 
index extremizes the total system damping. It has been shown in Ref. [El] 
that if one assumes that the disturbance is white noise entering the system 
at the controls with a magnitude N = R ^ and C = B then the results of both 
approaches yield the same result for the optimal calculation of the weight- 
ing matrix Q, however the optimal weighting matrix R is different. Using 
the same system dynamics as above, with Cov{n(t); nCx) = N6(t-T) unkno^m, 
one can summarize this approach as follows: 
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Cost function : ^ ^ / [/Ct)Qy(t)+u'^Ct)RuCt)]dt} 

uCt) ;l 

-T CEl-30) 


= Tr[CNC K] , (N: unknovm] 


If N = then 

o 

= kTr[BR‘^B'^K] 


C-El-31) 
CE1-32J 
CEl-35 ) 


where C : Ixn pseudo-inverse matrix of C,i.e., 


T -IT 

C=CC c) 


-1 T 


Optimal control law ; u(t) = -R B PCx(t3 


(El-34) 


where K: nxn unique positive definite solution of (El-21) 
KA + A^K - KBR'^b'^K + h'^QH = 0 (El-35) 

This optimal control law (22) can be computed without knowing N if 
Q and R are given. 


Covariance matrices (unkno'vn) : 

FX + XF^ + CNc"^ = 0 

where X = (x(t)x'^(t)}, F = A - BR“^B^K 

Y = {y(t)y^(t)}, EXh"^ 

U = (u(t)u^(t)} = E"^b'^KXKBR~^ 


(El-36) 
(El- 37) 

CEI-38) 


Note: These covariances are all unknown since N is unknom 


Auxiliary performance index : 

Total System Damping = TSD: the negative sum of the closed loop 

roots 

L* = Min Max {TSD} = Min Max {Tr [- (A-BR”^b'^K) ] } (El-39) 

Q R Q R 

or equivalently 

L* = Min Max {TrfBR ^b”^KJ} (’ . * A = constant matrix) (El-40) 

Q R 
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subject to KA + A^K - KBR'^b'^K + H^QH = 0 
Note: Cl) TSD is equivalent to the Relative Rate of Change of 

Error Sets (RERACES, see Ref. [1], Chapter IVj ^ and is 

generalized system response speed. 

C2) The extremum Crainimax) of TSD does not exist unless Q 

and R are appropriately constrained. 

Constraints on Q and R : [q[ = 1, |r| = CEl-41) 

where j(’)l : the determinant of C*) 

Optimal weighting matrices: 


Q = [hph'^ ! • chph’^)'^ 

R = . b'^KB ~ b'^KPKB 

m / - ' 

’'[b'^KB - B^KPKBl 

where FP + PF^ + BR~^B^ = 0 
F = A - BR'^b'^K 

Conclusions i 

Cl) When N is unknown, Q and R can be computed by the minimaximization 
of TSD with respect to Q and R under prescribed determinant con- 
straints . 

(2) The optimal Q matrix slows do\m the entire system response or 

error convergence speed by minimizing TSD. However, for an equi- 
valent white noise distrubance with a magnitude N inversely pro- 

© 

portional to the control weighting R, the minimization of TSD is 

equivalent to the minimization of J . CN jQjR). 

^ min e 


CEl-42) 

CE1-4S) 

CEl-44) 
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Note: a faster system response can be obtained only at the expense 

of the entire system accuracy, i.e., a larger TSD accompanfes 
a larger J . ( * ) • 

C3) The optimal R is chosen such that the system response speed (with 
a given Q selection) is maximized. 

(4) The optimal Q matrix equalizes the closed loop poles in fre- 
quency-wise, while the optimal R matrix tries to separate them. 
The trade-off between these two effects is controlled by p^. 

(5) is the only scalar parameter to be determined by designers. 

(6) The approach summarized here can be directly applied for the 
deterministic case in Summary 1 since Q and R can be determined 
independently of N and C matrices. 

Properties of the optimal weighting matrix determined by the OWEM design 
method (under determinant constraint), which can be of large significance 
in engineering applications, are the "Normalization and the Group-Wise 
Equalization". 

Normalization 

The state variables can be dimensionally different, and the weighted 
state and control variables are con^jared in the performance index. Hence, 
for a meaningful comparison, the weighted state and control variables should 
have the same dimensions. The resulting weighting matrices from the OWEM 
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design are inversely proportional to the output covariance matrix Y and the 
optimal control covariance. 

Group-Wise Equalization 

In the development o£ the OWEM design the followin^expression was 
obtained : 






= ks Y 


where 



C>0) 


Hence the Q is inversely proportional to the covariance matrix Y. 

The effect of the group-wise equalization constant k^ can be 
evaluated as follows: 




P/^P 


1 ^^ii ^ 

Q* = {qt.} = k„ Y*“-^ = k„ — 

11 ^ 3 Iy^i 


Yr^ : (i-j)th Cofactor of Y* = {yjyt} 


j-1 


qt . 

IJ 


ytyr 

11 


P 

L 

1=1 


Y*. 

-JLi- 

Iy*I 


ytyt 

11 


= k^; i=l,2. 




i = 1 : 

q*iiyf 

^ <^12 

V* V* + 

^2 a 

+ q? 

Ip 

y; 

II 

1 

CM 

\l 

H 

q^i y-fy-^ 

" ^22 

Y2>^1 

+ q* 

^2p 


= ^3 \ 

p-groups 

i = p : 


* V 

ypl ^ 

+ q* 
TP 

P P 




Similarly 


1=1 


11 


Cyty?) = k, E 


r 1 


* 

Y.. 

Ji 


1 = 1 iY"" 


ym 


= k. 


1 = 1 , 2 , 


This evaluation shows a group-wise equalization of the terms of 
{y*^Q y*}. 

In the initial approach all the constants are chosen to be equal. 

This corresponds to an Error Scaling Matrix (ESM) = 1. Making the ESM 
matrix not equal to the identity matrix corresponds to an error scaling 
factor (ESP) of the diagonal elements of the weighting matrix q^j . 

The elements of the ESM matrix are diagonal, hence, an easier inter- 
pretation is often possible. In general, the elements corresponding to 

1 

the quantities to be controlled determine the error weighting,while the 
elements corresponding to the derivatives of the quantities to be controlled 
influence the damping of that mode. In practical application, one should 
first use a ESM matrix with only elements that affect the quantity to be 
controlled. If mode damping increase is desirable, then elements of the 
ESM matrix that affect the mode damping can be introduced. This method, 
using a step by step approach, is preferable as the effects of the ESM matrix 
and Q matrix are not yet fully investigated. 
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APPENDIX E2 


THE OWEM METHOD APPLIED TO A SECOND ORDER SYSTEM 


A second order system equation is represented by the following 
differential equation 

xCt) = ax(t) + buCt} CE2-1) 

where 

xCt) = position 
x(t] = rate of position 
u[t) = control input 
a = rate damping 
b = control derivatives 
In the state variable form 


where 


x(t) = AxCt) + BuCt) 


xCt) « 


A = 


ct) 


x(t)' 

X£ (t) 


xCt)_ 

- 



r n 




0 1 


0 a 



u (t) = 6 

I 


CE2-2) 


(E2-3) 

CE2-4) 

(E2-5) 

(E2-6) 


230' 



An analytical approach using the example clarifies the method and 


yields a solution for the choices of the Q matrix. The closed loop 
system obtained by conventional optimal design method assigning values 
for qn, q 22 and r yields: 


xCt) 


+ v4^ + 2b / ^ t /ilL x(t)=0 (E2-7) 


"nie total system damping is (Ref. El-1) . 


TSD = 2b yaii- + b2 ai2 

r r 


(E2-8) 


The determinant constraint on weighting matrix Q becomes: 


iQj = 1 = qu q 22 


— (E2-9) 


This determinant constraint only specifies the lower limit of the mini- 
mized performance index value. Note that for decreasing qii, <122 and 
r, the value of the performance index decreases and the minimized value 
becomes meaningless. The determinant constraint is used to avoid this. 

Substituting Eq.(E2-3) into Eq. (E2-2) reduces the TSD to 


TSD = 


/a^-i- 


2b 





rqii 


(E2-10) 


231 




The minimization of TSD with respect to i elements leads to mini- 
mization of the performance index with determinant constraint. The mini- 
mization of TSD is obtained by differentiating (,E2-103with respect to qu. 
And TSD, qn and q 2 j arej 


TSD 


qn 


qzz 



CE2-U) 


(E2-12) 


(E2-13) 


Note that when the weighting factor Q is obtained, the OWEM method minimizes 
the performance index. 


J . 
min 



b3 

1 

tJ 


2 

XiCt) + 


1 

t5 

2 

b3 


2 

xz(t)+ r u^C^D) dt 

(E2-14) 
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Substituting the optimal weighting factors into fE2-l) gives the optimal 
closed loop system 



CE2-15) 


And the undamped natural frequency and the damping ratio C are re- 
spectively 






(E2-16) 


w = 
n 



(5 2-173 


The optimal control is given by substituting (E2-12) and (E2-15) into 
Cd 2-13)> 012-143, and from CD2-12 3. 


uCt3 = - ^ xCt3-C- f 

r y 




4a 

(br 3‘3 


-^kt) 


®2-18 3 


The solution of the Ricatti equation is as follows : 



CE2-193 
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Closed Loop System Properties as a Function of Control Weighting 


IThen the control weighting r goes to infinity, the closed loop system 
reduces the original basic one or its adjoint one 

x(t) + |a| x(t) = 0 (12-20) 

IVhen the control weighting r goes to zero, then 

10^ = oo (12-21) 

C = 0.866 (E2-22) 


which is more damped than that of the Butterworth poles for the second 
order system. 

Minimized Value of the Performance Index 
The minimized value of the performance index is: 


yJ . 

2 rain 


1 

y [kii ^10 2kx2 kai xio Xio + kzo ] 




+ 


2 

2r5 -Xio 


CE2-23) 
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when r = 0, 


It _ : 

2 rain 2 


CE2-24) 


Note that its value for conventional selection is given by 




which becomes zero when qaa = 0 • 
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APPENDIX FI 


TOE ROOT SQUARE LOCUS 


The Root Square Locus Method 
In the optimal control problem is minimized 


J = 


1 

2 



Ru) 


subject to 


x(tj = Ax(t) + Bu 
Y = Hx , H = I 

where x(t) : n dimensional state vector 

uCtj : ra diemsnioanl control vector (m^nj 
y(t) : p dimensional output vector (p<n) 

A,B and H: nxn, nxm and pxn constant matrices 

The solution of this problem yields the optimal control 

u* = -R“^B'^Kx(t) 

where K is a solution of the Riccati matrix equation 
KA + A^K - KBR'^B^K + Q = 0 
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CFl-1) 


CFl-2] 

CFl-3j 


CFl-43 


CFl-S) 



The optimal closed loop system is given by 


x(t) = Ax(t) + Bu* (t) 

= AxCt) + B(-R~^B'^K]xCt) (Fl-6) 

= (A-GK)xCt) CFl-7) 


where 


-1 T 

G = BR B 


(Fl-8) 


The characteristic equation is: 

|sl - A + BK| = 0 (Fl-9) 

The roots of the characteristic equation can be obtained by solving 
(Eq. Fl-93 with a priori -known K. The relationship between K,Q and R is 
given by C^q. Fl-5) , and is non-linear. 

The spectral form product of the characteristic Eq. Fl-9 

|sl - A + BK| 1-sI - (A-BK)'^I = 0 CFl-10) 

The characteristic equations of the optimal system and its adjoint 

X - Ax + = 0 (Fl-11) 

h'^QHx + X + a"^x = 0 


23 - 7 - 



-1 T 

IS - A BR B 

T T 

-H QH -IS -A 


CFl-1-2) 


and the following relationship is satisfied. 

jsi - A + BK[ I -si - CA-BK)"^! = IS-A BR'^B^ 

T T 

-H QH -Is -A 

(Fl-13) 

The roots can be obtained by solving (Eq. Fl-12} with respect to Q and R 
or graphically by making the "root square locus". Equation Fl-13 yields 
expressions relating the feedback gains K with A, B, H, and Q. 

Adding the control equation to (Fl-11) one obtains: 

X - Ax + BR'^b'^X = 0 

h"^QH + W a'^^X = 0 CFl-14) 

Ru + b”^X = 0 

-1 T 

and substituting u = -R B X into the top equation of (Fl-14), we can obtain 
the characteristic equation, 

T T 

-Is - A -H QH 0 

0 Is-A -B = 0 CFl-15] 

T 

B OR 


238 



In the transfer function form, this becomes: 


I + R 


-1 






= 0 


(Fl-16) 


where 


^(s5 = HCIs-Ar^ 





r (s) 

P 

Cs} 



ics) 

u 

m 


CFl-17) 


For single input system with R equal to a scalar p, the characteristic 
equation becomes : 




where 





CFi-is;) 


259 



The root square locus can be made by using similar rules as used for the 
root locus method, e.g. : 

N number of poles of 

p number of zeroes of 

Angle condition; 


(N-P3 = odd number 0 degree condition 
(N-p3 = even number 180° degree condition 


Asymptotes ; 


2mr 

2CN-P3 

N-p = odd 

n=0, ±1, ±2, 

2(N-p)-l 

2mr 

2CN-P) 

N-p = even 

n=0, ±1, ±2, 

2 (N-P) 


Center of asymptotes: always at origin 

2CN-P) poles go to infinity along asymptotes, as p^^O 
2 Cp) poles to zero, as p^^O 
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APPENDIX F2 


ROOT SQUARE LOCUS METHOD APPLIED TO A HELICOPTER CHOVER) 


The following is a summary from [Ref. FI ] , Consider a single 
input, dynamic system of the form: 

y(t)=AyCt) BuCt) CF2-1) 

where y(c) = n-dimensional state vector 

u(t) = m-dimensional control vector 
A, 3 = nxn and nxm constant matrices 
For an optimal control law, we seek to minimize the cost function. 

CO 

J(Q,R) = / [^*'^Ct3QyCt) + u’^(t)RuCt33 dt CF2-2X 

t 

o 

where Q,R = positive definite symme-cric constant 
weighting matrices. 


For this analysis, Q is chosen to be: 

qii 000 

0 ^,zz 0 0 

Q = 

0 0 qas 0 

0 0 0 q4L> 


and R = P a scalar 


(F2-3) 


CF2-4) 



The performance index (cost function) then becomes 


CO 

J = ~ f (yV + u’^’p^u) dt (F2-5) 

“ o 

The cnaracreristic equarion of the single input syscem is shorn. 
to be (Appendix FI) , 


I + p"^ C-s) (Sjj = 0 CF2-6) 

The characteristic equation then becomes; 


(-S), ^ (-S), ^ (-s)j 

"qn 

0 

0 

0' 

y, 

u 

Cs)~| 


0 

q22 

0 

n 

U. 

u 

(S) 


0 

0 

qa 3 0 

Zi. 

u 

Cs) 


0 

0 

d cfiti). 

_u 

(S) 


(F2-7) 


The change in roots of 'this equation^ resulting from varying a given 
parameter may be described by a root-square locus . The root-square 
locus technique is described in (Appendix FI) . 

A simplified longitudinal model for a helicopter is : 


(s-X^)u + g0 = Xg 

1 s 

-M^ u + s(s-M )9 = Mg Bj, 


(F2-8) 

(F2-9) 
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and may be expressed in state variable form as: 


V\ 


0 

1 

0 

0 


Vi 



0 

X 

u 

-g 

0 


Yi 

/3 


0 

0 

0 

1 


73 

74 


0 

M 

u 

0 

M 

q ^ 


71 


0 
0 

L“BjJ 


<u 


where 


_ *4 

7l 


’x' 

72 


X 

73 


0 

74_ 




(F2-10) 


(F2-11) 


From Eq» QF2-7) , the characteristic equation is 

1 * ^ [(qil ^ C-S) ^ (s) + qz2 ^ c-s) ^ (s: * q 3 3 ^ t-S) ^ Cs) 


, q*, ^ C-s) ^ Cs)] = 0 (F2-12) 

where ~ , — , — , -and ^ are given by; 

u ’ u ’ u ’ u * ^ 




CF2-13) 


CF2-14) 
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j 



Z3 - 0 - 


\s - X. * ^ U 




+M )s^ T X s + gM 
. q-' u q * t 


IF2-15) 


M_ „ J ^l5 \ 

s s-X +1 n Ig 

is uIm 3^J« 


s^-fX ■S'M }s^ ▼ X M s + ?M 1 

'• u q^ u q -' uj 


(F2-16) 


In order to determine the effect of q^^ on the roots of Equation F2-12, we 
apply the root-square locus technique. These cases will be examined^ 

First, weighting in the performance index will be applied only to yi, 
position. Next, the performance index will weight position and velocity. 
Weighting will then be applied to position, velocity and attitude. Then all 

four state variables will be weighted. 

Finally, is given a weighting for position and attitude. 

Case I 

Characteristic equation: 1 + ^ (-s) ^ (s) = 0 


J = 7 / Cqn yj +P dt 


(F2-17) 

CF2-18) 


The root-square locus is shown in Illustration F2-1. 
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Case II 


q y y q y y 

Characteristic equation: 1 + — (-s) (s) + — C-s) — (s) 

P P 


CF2-19) 


or: 


1 + yLSs) . ^ c-s} (i_s^ ^} = 0 

p2 u ■ U 


(F2-20) 


^ ~ jC +P^u^} dt 


(F2-21} 


The root-square locus is shown in Illustration F2-2. 


Case Ilia 

Characteristic equation: 

q y y q / y 

I ♦ -ii. ci Cs) C-s) t -f- ^s) ^C-s) 


= 0 


CF2-22) 


1 + ^ (s) ^ C-s} • 

P 

1 ♦ luLc.s^) !ls %S 

- ^ "b, \ 


11 


1= Cs^-M^s-j^g)(s^.M^s-3^g)J 


Is 


Is 


" J *^3 3y3 +P^u^ } dt 


CF2-23} 

CF2-24) 


The root square locus is shown in Illustration F2-B . 
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Case Illb 


Charact eristic equation: 


1 * V (-S) Cs3 . 


2 - U 


U 


2 U 


(.s) fs) * 2 m i (_ 3 J ( 3 ) 

P 


or 


1 + ii L_ • 

p 2 u Cs) u C-s) Cl + ^ C-s^3) 


1 + _33. {-_g2) 

qii ^ ^ 


\ 

«L * M^- V (— ’‘u * hT 

Is B, B, 


Is 


Is 


M 


B, 


M 


B, 


2 r^2 vf _ Is x;i „ 1^ /-I ^ 22 r „2 


K Cs'^-M^s- y q) Cs + M s- Y 

hs ®ls 




The root-square locus is shown in Illustrations F2-4. 


0 

(F2-2S} 


Case IV 

Characteristic equation: 

q y y 

1 + — ^ (-s) — ^ (s) 

p 



c-s) ^ Cs) 



— (-S) — Cs) 





Zi 

u 


C-s) 



0 


CF2-26) 


246 



or: 


q y y q 

1 + -LL -l_ (S) _J_ C-S}(;i + (-S^)) 

p2 u qii 


q 

qii 


CS-X + rr 

hs ^ 


q 

M ) (-S-X + M 3 (1 + C-s 3 
U-^ u M„ U-^ qa 3 


Is 


B 


Is 


M 


M. 


X2 (s^-M s - Cs^+ M s - Y 

IS =1 'I 


Is 


g3(l + c-s 

qii 


= 0 

(F2-273 


Th,e root square locus is shown in Illustrations {F2~53. 


\ 
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Conclusions 


1. In Case I, called position error control, where weighting is applied 
only to position in the performance criterion, the short period mode 
is not highly damped and the best possible phugoid damping is .707. 

2. In Case II, called velocity error control, where weighting is applied 
to both position and velocity, the damping of the attitude mode is 
still not high, but the position mode can be more damped than in Case I. 

3. In Case III, called acceleration error control, where weighting is 
applied to position, velocity and attitude, the damping of the attitude 
mode can be made greater than either of the two preceding cases. 

4. In Case IV, called accurate rate control, where weighting is applied to 
position, velocity, attitude, and attitude rate, the attitude mode is 
further damped. 

5. For precision flight path control with minimum (position) error, the 
velocity error control is minimum desirable. 

6. Acceleration and acceleration rate error control do not have a large 
effect upon the position error, but yield a reduction in control effort. 

7. In conclusion, for precision flight path control, only the short term 
dynamics for velocity, acceleration and accurate rate are important, hence 
one needs to consider only the "dynamic error control" for precision in 
flight path control. 
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Ilust^atiQU F2-1. Root-Square Locus - Weighting Applied to qn. Position 

^with varying — ^ 

P 

V ^ X 

tthe zeros are determined by the zeros of — or — ) 
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Illustration F2-2. Root Square Locus - Weighting 

Applied to qii (Position) and 
Ozz (Velocity) with varying qz 2 
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Illustration F2-3. Root Square Locus Weighting Applied to 
qH Cposition) and q^^ C^-'^titude) with Varying q^^ 
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Illustration F2-4. Root Square Locus Weighting Applied 

to qii (Position) qz 2 (Velocicy) 

* *q 

and qs 3 (Attitude) with varying — ^ 
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Illustration F2-5. Root Square Locus Weighting Applied 

to qii (Position), qz 2 (Velocity), 
qs3 (Attitude), and q 44 (Attitude Rate) 
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APPENDIX Gl 


REVIEW OF MODEL FOLLOWING BY LINEAR OPTIMAL CONTROL 
Introduction 
Given the system: 


x(t) = AxCt) + Bu(t3 


CGl-l) 


where: x(t) : n dimensional state vector, and output vector 

uft) : m dimensional control vector 
A, B : nxn, nxm and pxn constant matrices 
The "ideal" model is given by: 


yCt) = Ajj^ y(t) 

where : 


CGl-2] 


A is nxn constant matrix of the model 
m 

There are in principle two methods to provide a good match between the 
dynamic of the system and the dynamic of the "ideal model". Namely, 1. 
model in the performance index and 2. model in the system. 

Model in the Performance Index 

This approach is shown in Illustration Gl and uses the minimization 
of the norm: 


|x - xCt) I 1 


(Gl-3) 
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Illustration Gl-l. Model Following (Model in the Performance Index) 


Hence, the following performance index is defined. 


2V = Lim r CxCt) - A xCt))"^ Q Cx(t) - A xCt)) + u^Ct)Ru(t) dt 
0 


CGl-4) 


Using the Hamiltonian and maximum principle of Pontryagin (Ref. [Gl-2] , 
[Gl-9]), one obtains the optimal control law: 


u(t) = - (R'^b'^P + R ^B^QCA-Aj^) x(t) 


(Gl-S) 


where P is the nxn steady state solution to the Riccati equation: 


-PBR + a"^? + (A-A ) = 0 


(Gl-6) 
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and 


R = b'^QB+R 

A = A-BR"^B^QCA-A } 
m 

B = B 


CGl-7) 


When the dimensions o£ the model matrix are the same as that of system ma- 
trix and when the system and model are exactly matchable (Ref. G2 and G9)j 
the optimal control is of the following form: 

u*Ct) = -CA-Aj^) x(t) (G1-S-) 

Substituting Eq. (Gl-6) into [Gl-l) the optimal closed loop becomes: 

xCt) = A^ x(t) (Gl-9) 

which has the dynamic of the model. This exact matching is achieved (in 
general) as IqI/IrI becomes large. (Ref. [Gl-2]). 

IVhen the model and system are exactly matchable, there is no need in 
using optimal control and the optimal control *(t) can be derived from 
algebraic considerations (Ref. [Gl-6]). The algebraic conditions for "perfect" 
following are derived in Ref. [Gl-6] . In order that the dynamics of the model 
and system be the same, it is reqiiired that: 

x(t) = A^ x(t) (Gl-10) 
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Substituting CGl-8) into Gl-l), 


x(t) = Ax(t) + Bu(t) (Gl-11) 

hence 

Bu(t) = x(t) CGl-12) 

Multiplying both sides of (Gl-10) by B Cwhere B is the pseudoinverse of B) 
one obtains 

u(t) = b'*’{;a^-A] xCt] (Gl-13) 

It has been shown in Ref. [G1-6J that necessary and sufficient conditions 
for perfect matches between the dynamics of the plant and that of the model 
are: 

[B b"‘-I][A^-A] = 0 CGl-14) 

This equation indicates that for u(t) (given by Gl-13) to exist, there must 
be a great structural similarity between the system matrix A and the model 
matrix L. Moreover, B must have sufficient rank. However, it is clear that 
for this limited class of problems, this approach involves much less computa- 
tion than the quadratic cost integral approach. 

Model in the System 

This approach is shown in Illustration G2 and is based on the minimization 
of the square error between the output of the given plant and the output of the 
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desired plant model; thus forcing the output of the given system to 
follow the model. 



Illustration Gl-2. Model Following (Model in the System) 


Given the system; 


XgCt) = A^x^(t) + B^u^Ct) 


(Gl-15) 


where : 

x^(t) : n dimensional state vector, and output vector 
u^(t) ; m dimensional control vector (m<n) 

A^,B^ : nxn, nxm and pxn constant matrices 

The model is given by: 


2SS 



where : 


x^(t) ; n dimensional state vector, and output vector 
u^Ct) : m dimensional control vector (nKn) 

A , B : nxn, nxm and pxn constant matrices 

In general, the dimension o£ the system's output may be different from that 
of the model (Ref. [Gl-10], Gl-7] and [Gl-2]). The appropriate performance 
index is defined as : 

T 

^ ^ I ~ ^ ~ ^ u(t)]dt 

(Gl-17) 

where : 

Q ; nxn positive-semi definite matrix 
R : vxv positive-definite matrix 

The performance index J is to be minimized with respect to the control 
uft). This problem can be treated as tracking problems or servo -problems. 
CRef.[Gl-8], Gl-3] , [Gl-10]). There are basically two different approaches 
to the problem. One is established by Ref. [Gl-6]; in that work algebraic 
conditions for perfect following are derived. The second is developed by 
Ref. [Gl-1], and[Gl-9] and others. The development in[ Ref . Gl-l] is basically 
the same as in jRef . Gl-9] with the exception that the input to the model is a 
zero input. 

^Vhen given the system CGl-15]) and the model CGl-16), and suppose that u^^ 
can be produced by ordinary differential equations such that; 

u = Du (Gl-18) 

mm ^ ^ 
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then one can define the enlarged state vector x(t): 


x(t) = 


X (t) 

P 


X Ct) 

ra 


u (t) 


The dynamic of the new system is given by: 


where : 


where : 


x(t) + Bu(t) 


A' 0 0 

s 


A = 


m 


m 


B 


B, 

0 

0 


A is a C2n+m) x(2n+m) matrix 
B is a (2n+m)x v matrix 


CGI -193 
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Considering the new system a regulator problem (Ref. [Gl-1] , [Gl-10] ) , 
one can define the following performance index: 

T ^ 

J = ^ [x'^Ct) Q x(t) + u’^(t) R u(t)] dt (Gl-21) 

where R is a constant positive definite matrix, and J is to be minimized 
with respect to u(t). In order that J is meaning full, Q must be defined 
in appropriate ways. Since the error between the system and the model is 
to be minimized, the weighting matrix Q is defined as follows: 

" Q -Q 0“ 

Q = -Q Q 0 (Gl-22) 

0 0 0 _ 

It is easy to verify that with the definition Eq. (Gl-22) substituted into 
Eq. (Gl-21), the performance index in Eq. (Gl-3) is identical to that in Eq. 
(Gl-9). Hence, the model following is formulated by defining the performance 
index. Using the Hamiltonian for the system, Eq. (Gl-6) and the performance 
index Eq. (Gl-9) and applying Pontryagin’s maximum principle, one obtains the 
following optimal control (Ref . [Gl-1] , [Gl-7] and others). 

u*(t) = * Pi3%Ct)] 

(Gl-23) 
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or 


u*(t) = -K3X3C-t) + 


where 


-1 T 

K = R B P,, 
s s 11 

-1 T 

K = -R ■'b 'p,_ 
m s 12 


-1 T 

K = -R B P,, 
V s 13 
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P^^ is the solution to ricatti equation, and is nxn syinraetric matrix. 

^11 ^lA ^ ^s^^l ^ Q ' ^ 11^^11 = ° 

The solution P^^ o£ Eq. (^1-26) is independent of the other two gain- 
matrices : Pj^2> ^23' matrices (Eq. Gl-25) are the solutions to 

the following algebraic equations : 

'■l2 ♦ ^2 * ■’lA = Q (=1-27) 

■’is * fl3 * ’’ij” ■ 

where 

-1 T 

A = B R B 
s s 
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The solution of Equation (G1-2S) depends on D which is undesirable. 
Therefore, it usually makes D equal to zero which means designing the system 
to step input (see for example. Ref. [Gl-9]). 

For linear systems with time-'invariant parameters, it is very desirable 
to find a solution (P^^ matrices) which is time-invariant which means constant 
gams. The question which arises is if there exists the limit of the perfor- 
mance index J as T approaches infinity. It is well known that sufficient con- 
dition for this limit to exist is that the system (Eq. Gl-1) is completely con- 
trollable. Unfortunately the system (Eq, Gl-1) is uncontrollable by defini- 
tion, so that the existence of the limit of (Eq. Gl-3) is not guaranteed. 

In Ref. Gl-5 is shorn that partially controlled systems can be divided 
into two parts, the part being independent of each other. Hence, the Ricatti 
equation has the well known steady state solution. Similarly as in the de- 
velopment of the model in the perfo3rmance index in Ref. Gl-1 and Gl-3. 

Conditions are established for perfect model following. In Ref. Gl-3 
is taken the difference between the system and the model equations (Gl-15, 
Gl-16) and substituted to that the optimal control given by Eq. (Gl-25), we 
have 


X (t) - X (t)) = (A - B K )x (t) + (B K -A )x (t) 
s m s ss s ^nim m m 


+ (B K - B )u (t) 


(Gl-29) 


263 



If one equates: 


A - BK =-BK + A 
s s s pm m 


B K = 
s V m 


Then from Eq. (Gl-29] 


CGl-30) 


(x^Ct) - x^Ct3) = (Ag-B^K^Kx^Ct) - (Gl-31) 

and may be chosen to make the system less sensitive to variations of 
parameters. Up to this point, this condition (Eq. Gl-30) for perfect following 
seems to be equivalent in both approaches. Conditions require that is 
invertible, which often is not the case. 

In the case where is not invertible. Ref. [Gl-7] canonical forms 
(for derivation, see Ref. Gl-12) in order to solve (Eq. Gl-30). 

It has been said that for any pair (A^, B^) there exists a transformation 
T such that: 


-1 T 
T Ag = 


^ J f 

- I- - 

A / B 


} n-m 
} m 



B 


0 } n-m 

T } m 
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where 0 is null matrix. However, such transformation does not exist in 
general (see Ref. Gl-13) . Hence, Eq. (Gl-30) does not have a solution for 
every system and model, and this reduces the number of cases which can be 
solved. 

In Ref. (Gl-3) is derived the algebraic solution for model following; 
following that derivation for model following a command input and using the 
previous notation, one obtains : 

A 

X (tl = A x(t) + B u (t) 

A 

Where is command input and B as a constant command control matrix. 

For perfect matching: 

X (t) = X (t) 

and also its first derivatives should be equal. This yields the control law: 
u(t) = (B^KAj^-Ap x(t) r B^'b u^(t) 

The first derivatives of and x^(t) will be equal for all x(t) and 

if 

[BB^-I] [A^-A^] x(t) = [BB^-I]B u^(t) 
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which yields the conditions 


[B B - I] CA -A ) = 0 
*■ s s ^ m 

. - (Gl-32) 

[B B - I] B = 0 
'■ s s 
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